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ABSTRACT 

With the goal of investigating the degree to which the mid-infrared emis- 
sion traces the star formation rate (SFR), we analyze Spitzer 8 /xm and 24 /xm 
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data of star-forming regions in a sample of 33 nearby galaxies with available 
HST/NICMOS images in the Paa (A1.8756 fim) emission line. The galaxies 
are drawn from the SINGS sample, and cover a range of morphologies and a 
factor ~10 in oxygen abundance. Published data on local low-metallicity star- 
burst galaxies and Luminous Infrared Galaxies are also included in the analysis. 
Both the stellar-continuum-subtracted 8 /zm emission and the 24 /xm emission 
correlate with the extinction-corrected Fact line emission, although neither re- 
lationship is linear. Simple models of stellar populations and dust extinction 
and emission are able to reproduce the observed non-linear trend of the 24 fim 
emission versus number of ionizing photons, including the modest deficiency of 
24 /im emission in the low metallicity regions, which results from a combination 
of decreasing dust opacity and dust temperature at low luminosities. Conversely, 
the trend of the 8 /xm emission as a function of the number of ionizing photons 
is not well reproduced by the same models. The 8 /im emission is contributed, in 
larger measure than the 24 /zm emission, by dust heated by non-ionizing stellar 
populations, in addition to the ionizing ones, in agreement with previous find- 
ings. Two SFR calibrations, one using the 24 /im emission and the other using 



a combination of the 24 /xm and Ha luminosities (iKennicutt et al.l l2007al ). are 
presented. No calibration is presented for the 8 /im emission, because of its 
significant dependence on both metallicity and environment. The calibrations 
presented here should be directly applicable to systems dominated by on-going 
star formation. 

Subject headings: galaxies: starburst - galaxies: interactions - galaxies: ISM - 
ISM: structure 



1. Introduction 

The mult i- wavelength galaxy surveys of unprecedented angular resolution recently made 
available by combined space (HST, Spitzer) and ground-based observations are providing 
for the first time the tools to cross-calibrate star formation rate (SFR) indicators at different 
wavelengths, and to test the physical assumptions underlying each indicator. 

Easy accessibility has traditionally favored the use of the ultraviolet (UV) stellar contin- 
uum and of the optical nebular recombination lines as SFR indicators, the former mainly in 
the intermediate-high redshift regime (as it gets redshifted into the optical observer frame) 
and the latter mostly in low-redshift surveys. Both indicators only probe the stellar light that 
emerges from a galaxy unahsorhed by dust. The UV is heavily affected by dust attenuation. 
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and numerous efforts have attempte d to find general tools to mitigate the effects of dust on 



this important SFR indicator (e.g., 


Calze 


tti. Kinnev & Storchi-Beremann 


1994 




Kennicutt 


1998a 


Meurer. . 


ieckman & Calzetti 


1999; 


Hopkins et al. 


2001; 


Sullivan et al 


. 2001 


Buat et al. 


2002 




2OO5I; 


Bell 


I2OO3; 


Hopkins 


2004 


: Salim et al. 


2007 


). Cross-calibrations with optical re- 



combination lines and other indicators have also attempted to account for the ~10 times or 
more longer stellar t imescales probed by the UV relative t o tracers of ionizing photons (e.g., 
Sullivan et al.ll200ll ; iKong et al.ll2004l : ICalzetti et al.ll2005l ). Among the recombination lines. 
Ha is the most widely used, due to a combination of its intensity and a lower sensitivity 
to dust attenuation than bluer nebular lines. Although to a much lesser degree than the 
UV, the Ha line is still affected by dust attenuation, plus is impacted by assumptions on 



the underlying 


stellar absorption and on the form of 


;he high end of t 


le ste 


iar initial mass 


function (e.g.. 


Calze 


tti. Kinnev & Store 


li-Bergmann 


1994; 


KennicuttI 


1998a 


Hopkins et al. 


2001; 


Sullivan et al. 


2001; 


Kewlev et al. 


2002; 


Rosa-Gonzalez. Terlevich & Terlevich 


2002 


)■ 



Infrared SFR indicators are complementary to UV-optical indicators, because they 
measure star formation via the dust-absorbed stellar light that emerges beyond a few 



the IRAS times (e.g.. 


Lonsdale Persson & Helou 


1987; 


Sauvage & Thuan 


1992!). interest in the this wavelength 



cent times b y the discovery of subraillimeter-emitting galaxy populations at cosn i ological dis- 



tances (e.g.. ISmail. lyison fc Blainlll997l ; iHughes et al.lll998l ; iBarger et al.lll998l ; lEales et al. 



I999I ; IChapman et al.ll2005l ). In dusty starburst galaxies, the bolometric i nfrared luminosit y 



LjR in the ~3-1100 /xm window is directly proportional to the SFR ( e. g. . lKennicuttlll998al ) . 
However, even assuming that most of the luminous energy produced by recently formed stars 
is re-processed by dust in the infrared, at least two issues make the use of this SFR indicator 
problematic. (1) Evolved, i.e., non-star-forming, stellar populations also heat the dust that 
emits in the IR wavelength region, thus affecting the calibration of SFR(IR) in a stellar 



popula tion-dependent manner (e.g.. iLonsdale Persson fc Heloulll987l : lHeloulll986l ; iKennicutt 



1998al ). (2) In intermediate/high redshift studies, the bolometric infrared luminosity is often 
extrapolated from measurements at sp arsely sampled wavelengths, most often in the sub - 



mm and radio observer's frame (e.g., ISmail. Ivison fc BlainI 119971 ; IChapman et al.l l2005l ) 



and such extrapolations are subject to many uncertainties. 

The interest in calibrating monochromatic mid-infrared SFR diagnostics stems from 
their potential application to both the local Universe and intermediate and high redshift 



galaxies observ ed with Spitzer and future infrared/submillimeter missions (jPaddi et al.ll2005 



Wu et al.ll2005l ). One such application is the i nvestigation of the scaling laws of star form ation 
in the dusty environments of galaxy centers (jKennicuttlll998bl ; IKennicutt et al.ll2007al ). The 
use of monochromatic (i.e., one band or wavelength) infrared emission for measuring SFRs 
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offers one definite advantage over the bolometric infrared luminosity: it removes tlie need 
for liiglily uncertain extrapolations of the dust spectral energy distribution across the full 
wavelength range. Over the last few years, a number of efforts have gone into investigating 
the potential use of monochromatic infrared emission for measuring SFRs. 

Early studies employing ISO data have not resolved whether the warm dust and aromatic 
bands emission around 8 fim can be effectivel y used as a SFR indi cator , since different conclu- 
sions h ave been reached by different authors. iRoussel et al.l (120011 ) and lForster Schreiber et al. 
(120041 ) have shown that the emission in the 6.75 /zm ISO band correlates with the number 
of ioniz ing photons (SFR) in galaxy disk s and in the nuclear regions of galaxies. Con- 
versely, iBoselli. Lequeux &: Gavazzil (120041 ) have found that the mid-IR emission in a more 
diverse sample of galaxies (types Sa through Im-BCDs) correlates more closely with trac- 
ers of evolved stellar populat ions not linked to the current star formation. Additionally, 
Haas. Klaas fc Bianchil (12002h find that the ISO 7.7 /im emission is correlated with the 
850 /im emission from galaxies, suggesting a close relation between the ISO band emis- 
sion and the cold dust heated by the general (non-star-forming) stellar population. This 
divergence of results highlights the multiplicity of sources fo r the emission at 8 /xm (e.g., 
Peeters. Spoon fc Tielend l2004l : iTacconi-Garman et al.ll2005l ). as well as the limits in the 
ISO angular resolution and sensitivity for probing a sufficiently wide range of galactic con- 
ditions. 



The emission in the 8 /i m and other MIR bands is generally attributed to Polycyclic Aro- 
matic Hydrocarbons (PAH, Leger fc Pugetlll984 : Sellgren 1984 : AUamandola. Tielens &: Barker 
19851 : ISelkren. Luan fc Wernei]ll99oh 7 large molecules transiently heated by single UV and 



optical photons in the general radiation field of galaxies or near B stars ( Li fc Draine 



Haas. Klaas fc Bianchi 



Wu et al. 



2005 



Houck et al. 2004 



2002 



2002l:lBoseUi. Leaueux fc Gavazzl20oi : IPeeters. Spoon fc TielenJ2004j : 



Mattioda et al. 1 120051). and which can be destroyed, fragmented, or ion- 



ized by harsh UV photon fields (IBoulanger et al.l Il988l . Il990l : iHelou. Ryter fc Soifed Il991 



Pety et al. 20051 ) . Spitzer data of the nearby galaxies NGC300 and 



NGC4631 show that 8 /im emission highlights the rims of HII regions and is depressed 
inside the regions, indicating that the PA H dust is heated in the PDRs surrou nding HII 



regions and is destroyed within the regions (IHelou et al. 



2004 



Bendo et al.ll2006l ). Analysis 



of the mid-IR emission from the First Look Survey (iFang et al.ll2004j ) galaxies shows that 
the correlation between t he Spitzer 8 /trn . band emission and tracers of the ionizing photons 
is shallower than unity (IWu et al.l l2005l ). in agreement with the correlation s observed for 



HII r egions in the nearby, metal-rich, star-forming galaxy NGC5194 (M51a ICalzetti et al. 
2005h . 



The 24 /im emission is a close tracer of SFR in the dusty center of NGC5194 (ICalzetti et al, 
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2005h and in NGC3031 JPerez-Gonzalez et all [2OO6I 1 The general applicability of this 
monochromatic indic ator has so far been explored only for a sra all number of cases, mostly 
bright galaxies (e.g.. lWu et al.ll2005l : lAlonso-Herrero et al.ll2006l ). A potential complication 
is that most of the energy from dust emerges at wavelengths longer than ~40-50 fim (see 



Dale et al.ll2006l . and references therein). Thus the mid-IR does not trace the bulk of the 
dust emission, and, because it lies on the Wien side of the blackbody spectrum, could be 
sensitive to the dust temperature rather than linearly correlating with source luminosity. 

This study investigates the use of the Spitzer IRAC 8 fim and MIPS 24 /im monochro- 
matic luminosities as SFR indicators for star forming regions in a subsample of the SINGS 
galaxies (SINGS, or the Spitzer Infrare d Nearby Galaxies Survey, is one of the Spitzer 
Legacy Programs, iKennicutt et al.ll2003l ). Star-forming regions in galaxies represent a first 
stepping-stone for characterizing SFR indicators, as they can be considered simpler entities 
than entire galaxies. 

We also extend our analysis to include both new and published integr ated (galaxy-wide) 
data on local low-metal licity starburst galaxies (lEngelbracht et al.l 120051 ) and Luminous In- 
frared Galaxies (LIRGs. lAlonso-Herrero et al.ll2006l ). These data are used to explore whether 
the relationships derived for the star-forming regions that constitute our main sample are 
applicable to starburst-dominated galaxies as a whole. A future paper will investigate the 
viability of the mid-infrared lum inosities as SFR tracers for more general classes of galaxies 
kennicutt fc Moustakaslbood ). 



The Spitzer observations are coupled with near-infrared HST/NICMOS observations 
centered on the Paschen-a hydrogen emission line (Paa, at 1.8756 /im), and with ground- 
based Ha observations obtained by the SINGS project. The hydrogen emission lines trace the 
number of ionizing photons, and the Paa line is only modestly impacted by dust extinction. 
Furthermore, the Paa and Ha lines are sufficiently separated in wavelength that reliable 
extinction corrections can be measured (iQuillen fc Yukital 1200 ll ). Because of its relative 
insensitivity to dust extinction (less than a factor of 2 correction for the typical extinction 
in our galaxies. Ay ^5 mag), Paa represents a near ly unbiased tracer of the current SFR 
over a timescale of about 10 Myr (IKennicuttl Il998al ) . The access to Paa images to use as 
a yardstick for calibrating the mid-infrared emission is the basic motivation for the present 
work. 

The present paper is organized as follows: Section 2 introduces the sample of local star- 
forming galaxies from SINGS; Section 3 presents the data, while the measurements used 
in the analysis are pre sented in Section 4 . Sect ion 5 briefly introduc es the low metallicity 



starbu rst galaxies from lEngelbracht et al.l (120051 ) and the LIRGs from lAlonso-Herrero et al. 
( 120061 ). The main flndings are reported in Section 6, and the comparison with models is 
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made in Section 7. Discussion and a summary are given in Sections 8 and 9, respectively. 
Details on the models of dust absorption and emission are in the Appendix. 



2. Main Sample Description 



The SINGS sample of 75 galaxies (IKennicutt et al.ll2003l ) was used as our baseline sample 
for which HST observations in the infrared were either obtained as part of our project or 
retrieved from the HST archive (see details in section 3.2). The only criterion required for 
a SINGS galaxy to be observed with the HST was to have a redshifted Paa emission within 
the transmission curve of one of the NICMOS narrowband filters. A total of 39 galaxies, 
or 52% of the SINGS sample, were observed in the Paa line (example in Figured]). The 
HST/NICMOS-observed galaxies are fully representative of the SINGS sample as a whole, 
in terms of morphological types, range of metallicity, and SFRs. 

The infrared data of 4 of the 39 galaxies show non-recoverable problems (see section 3.2 
for additional explanation); two more galaxies, MSlDwA and DD0154 do not show either 
optical line emission or mid-IR dust emission in the region imaged in the near-infrared with 
HST. All six galaxies were discarded from the current analysis, thus leaving a net sample of 
33 galaxies. Table [H lists the main characteristics of the 39 galaxies, separating the discarded 
ones from the remainder of the sample. 

The 33 galaxies are divided in three groups according to their oxygen abundance: high 
metalhcity galaxies (12+log(0/H)>8.35), medium metalhcity galaxies (8.00<12+log(O/H)<8.35), 
and low metallicity galaxies (12+log(O/H)<8.00). The two sets of disk- averaged oxygen 



abundance values listed in Table [J differ systeni atically by about 0.6 dex (IMoustakas et al. 
20071 ) . As described by IMoustakas et al.l (120071 ) , the set of lower numb ers for the oxygen 



abundance is roughly tied to the electron temperature abundance scale (jPilyugin &: Thuan 
20051 ) . whi le the higher abundance set is based on stellar popu lations plus photoionization 
modelling (jKobulnicky &: Kewley 2004; Kewley fc Dopita 2002). The difference between the 
two scales is due to a as-yet unidentified systematic zeropoint offset, and the 'true' oxygen 
abundance should lie somewhere between the two listed values; however, the relative ranking 
of abundances on either of the scales should be fairly accurate. On this basis, we assign a 
galaxy into a metallicity bin based on the average of the two values. Metallicity gradients 
across galaxies are likely of little impact in our analysis. The observations probe the inner 
~0. 8-5.1 kpc, depending on the distance; typica l metallicity variations over these region 
sizes are less than ~0.3 dex for our spiral galaxies (IMoustakas et al.ll2007l ). and therefore are 
not expected to play a significant role in our results. 
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Within the area imaged by the HST/NICMOS for each galaxy in the main sample 
(Table [1]), regions of star formation are identified and their fluxes measured over typical 
sizes of ~200-600 pc (section 4 and Figured]). These regions are termed here 'HIT knots', 
and they are far simpler units, in terms of stellar population and star formation history, than 
whole galaxies. The HII knots in this study cannot be considered individual HII regions in the 
strict meaning of the term. Limitations in angular resolution, as discussed in section 4, force 
us to consider areas within galaxies which may be populated by multiple HII regions. The 
main requirement is for such areas to be local peaks of current star formation, as determined 
from hydrogen line or infrared emission. The ionizing populations in these regions can be 
approximated as having comparable ages, and more evolved stellar populations do not tend 
to dominate the radiation output. Although caution should be used when deriving a star 
formation rate for quasi-single-age populations, the investigation of simpler, star-formation- 
dominated structures should offer better insights than whole galaxies on the strengths and 
weaknesses of the mid-infrared SFR indicators of interest here. 



Details on t he low-metallicity sta r burst galaxies from lEngelbracht et al.l (120051 ) and on 
the LIRGs from lAlonso-Herrero et al.l (l2006l ) are given in section 5. 



3. Observations and Data Reduction 

3.1. Spitzer IRAC and MIPS Imaging Data 

Spitzer images for the galaxies in Table [H were obtained with both IRAC (3.6, 4.5, 5.8, 
and 8.0 /im) and MIPS (24, 70, and 160 /xm), as part of the SINGS Legacy project, between 
~March 2004 a nd ^August 2005. A d escription of this project and the observing strategy 
can be found in iKennicutt eTaP J2003h . 



Each galaxy was observed twice in each of the four IRAC bands, with a grid covering 
the entire galaxy and the surrounding sky. The observing strategy allowed a separation of a 
few days between the two observations to enable recognition and exclusion of asteroids and 
detector artifacts. Total exposure times in each filter are 240 s in the center of the field, 
and 120 s along the grids' edges. The SINGS IRAC pipeline was used to create the final 
mosaics, which exploits the su b-pixel dithering to better sample the emission, and resamples 



each mosaic into 0.75" pixels ( iRegan et al.ll2004j ). The measured 8 /im PSF FWHM is, on 



average, 1.9", and the 1 a sensitivity limit in the central portion of the 8 fim mosaic is 
1.2x10"^ Jy arcsec"^. 

As the interest in this paper is in using the dust emission at mid-infrared wavelengths 
(8 /im and 24 /im) as SFR tracers, we need to remove the stellar continuum contribution 
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from the 8 /x m images. This co ntribution is, in general, small in high metallicity, dusty 



galaxies (e.g.. lCalzetti et al.ll2005l ). but can become significant in lower metallicity, and more 



dust -poor galaxies. 'Dust-emis sion' images at 8 u m are obtained by subtracting the stellar 
contribution using the recipe of iHelou et al.l (120041 ): 



dust{^) — fsfimi^) — P fsSfimi^) , (1) 

where the coefficient (3 is in the range 0.22-0.29, as determined from isolated stars in the 
galaxies' fields. Visual inspection of the stellar-continuum subtracted images suggests that 
this approach is fairly accurate in removing stellar emission; occasional foreground stars 
located along the galaxies' lines of sight are in general removed by this technique. Although 
the 3.6 fim images can include, in addition to photospheric emission from stars, a component 
of hot dust emission, this component is unlikel y to have an inapact beyond a few percent on 
the photometry of the dust-only 8 fim images (jCalzetti et al.l 120051 ). 



MIPS observations of the galaxies were obtained as scan maps, with enough coverage 
to include surrounding b ackground in additio n to t he galaxy. The red uction steps for MIPS 
mosaics are described in iGordon et al.l (120051 ) and iBendo et al.l (120061 ). At 24 fim, the PSF 
FWHM is ~5.7", and the 1 a detection limit is 1.1x10"^ Jy arcsec"^. The MIPS images are 
considered 'dust' images for all purposes, as contributions from the photospheric emission of 
stars and from nebular emission are negligible (a few percent) at these wavelengths. 



3.2. HST Imaging Data 

The main advantage of using near-infrared narrowband imaging, rather than spec- 
troscopy, is the potential of capturing, in principle, all of the light in the Pact line, thus 
enabling a more secure measurement of the total line emission from the targets. The 
HST/NICMOS narrowband filters of interest here have ~1% band-passes, that can eas- 
ily accommodate gas line emission with a few hundred km/s shift relative to the galaxy's 
systemic velocity. 

Most of the HST/NICMOS observations for the galaxies in our sample come from the 
HST SNAP program 9360 (P.I.: Kennicutt). For 9 of the galaxies, archival HST data were 
used, from programs GO-7237 and SNAP-7919. 

Observations for SNAP-9360 were obtained with the NIC3 camera, in the narrowband 
filters F187N, F190N (Paa emission line at restframe wavelength A=1.8756 /xm and adjacent 
stellar continuum), and the broadband filter F160W. The NIC3 camera has a field of view of 
51", and observations were obtained with 4 dithered pointings along a square pattern with 
0.9" sides, to better remove cosmic rays and bad pixels. Thus, NICMOS observations imaged 



- 9 - 



the central <1 arcmin of each galaxy. The NIC3 0".2 pixels undersample the NICMOS PSF, 
although this is not a concern for the diffuse ionized gas emission. On-target total exposure 
times were 640 s, 768 s, and 96 s, for F187N, F190N, and F160W, respectively. 

The data were reduced with the STScI IRAF/STSDAS pipeline calnica, which removes 
instrumental effects, bad pixels, and cosmic rays, and produces images in count-rate units. 
The removal of the quadrant-dependent 'pedestal' was done with the IRAF/STSDAS routine 
pedsub. The four dithered exposures were combined with the IRAF/STSDAS mosaicing 
pipeline calnicb. 

For our analysis, only the two narrowband images are used, and the emission line-only 
images are obtained by subtracting the continuum-only images, rescaled by the ratio of the 
filters' efficiencies, from the line+continuum image. Program 9360 was executed after the 
NICMOS Cryocooler System (NCS) had been installed on the HST, providing a detector 
quantum efficiency about 30% higher in the H-band than during pre-NCS (i.e., pre-2002) 
operation^. This is an important difference when comparing depths of SNAP-9360 with 
those of the archival NICMOS images, which were obtained pre-NCS. The average 1 a 
sensitivity limit of the continuum-subtracted image is 6.4x10"^'' erg s~^ cm~^ arcsec"^. In 
units that will be easier to relate to the analysis performed in this paper, our 1 a limit for 
a specific Paa luminosity measured in a 13"-diameter aperture is 2.83x10^^ erg s~^ kpc^^; 
in a 50"-diameter aperture, the 1 a limit is 1.04x10^^ erg s~^ kpc~^ . 



The archival NICMOS data from HST snapshot program 7919 are described in lBoker et al 



(119991 ). Here we summarize the main differences with SNAP-9360. Data for the SNAP- 79 19 
were obtained with a single pointing (and a single integration) of the galaxy's center with 
the NIC3 camera. One narrowband filter (F187N or F190N depending on redshift) and the 
broadband F160W filter were used, for ~768 s and 192 s, respectively. We re-processed 
the archival images through calnica, to improve the removal of instrumental effects and of 
co smic rays by using a more recent version of the calibration pipeline than the one used 



Boker et al.l (119991) : the quadrant-dependent pedestal was removed with pedsub. As in 

Boker et al. rjl999h . the rescaled broadband filter is used for removal of the underlying stellar 



continuum from the image containing the Paa emission line. The images from SNAP-7919 
are deeper than in SNAP-9360, with an average 1 a sensitivity limit of the continuum- 
subtracted image of 3.5x10^^'' erg s~^ cm~^ arcsec"^. 

Broadband filters may not provide the optimal underlying stellar continuum signature, 
especially if uneven dust extinction in the galaxy produces color variations within the fil- 



^The Near Infrared Camera and Multi-Object Spectrometer Instrument Handbook, version 9.0, E. Barker 
et al. eds., 2006, STScI 



- 10 - 



ter's bandpass. To check the impact of this potential effect, we have compared observa- 
tions of galaxies in common between the SNAP-9360 and SNAP-7919 programs: NGC3184, 
NGC4826, NGC5055, and NGC6946 (images of NGC0925 are also present in both programs, 
but the pointings are only partially overlapping, and are sufficiently different that both im- 
ages are used in our analysis, see Tabled]). For SNAP-9360, two narrowband images are 
available, thus yielding a 'cleaner' line image. Comparison of continuum-subtracted images 
in both programs for regions in the common galaxies yields differences in the Paa photom- 
etry in the range 10%-30%, which is in general well within our random uncertainty for the 
Paa measurements (section 4.2). 



The NI CMOS archiva l data for NGC5194 (HST program 7237) are described in lScoville et al. 



(I2OOII ) and ICalzetti et al.l (120051 ). The main difference with the data in 9360 is that the 
NGC5194's image is a 3x3 NIC3 mosaic that spans the central 144" arcsec^. Each pointing 
was observed in both F187N and F190N, with 128 s exposure times. The sensitivity is vari- 
able, being lower at the seams of the 9 images that form the mosaic. The average 1 a sensitiv- 
ity limit of the continuum-subtracted image for this galaxy is 1.8x10"^^ erg s~^ cm~^ arcsec"^. 

The HST/NICMOS observations are the shallowest in our sample when compared to the 
other images, and represent the true limitation to our analysis. On the other hand, Paa mea- 
surements offer an opportunity to obtain a nearly unbiased measure of the number of ionizing 
photons produced in a region, as it is only weakly affected by dust extinction. An extinction 
as large as 5 mag at V produces an extinction of 0.73 magnitudes at Paa, i.e., roughly a 
change of a factor of 2 in the line intensity (Figure [2]), for foreground screen dust geometry. 
Still, we combine the Pact measurements with complementary measurements at Ha to cor- 
rect the line emission for the effects of dust. We adopt a metallicity-dependent intrinsic ratio 
HQ;/Paa=7.82, 8.45, and 8.73 for the high, medium, and low metallicity subsamples, respec- 
tively, which correspond to elec tron temperatures Tp=7,000 K, 10,000 K, and 12,500 K for 



the H II knots (for ne=100 cm lOsterbrock fc Ferlandll2006l : iGarnett. Kennicutt fc Bresolin 
2004). We also adopt an extinction curv^ with differential value k(Ha;)— k(PaQ;)=2.08 



flFitzpatrick Ill986l : iLandini et al.lll9841 ). 

Four of the galaxies discarded from our sample (Table [T]) present an array of problems 
mainly in their NICMOS observations. The F187N image of NGC0024 is heavily affected by 
cosmic ray persistence, which has caused the effective noise level of the frame to be about 
7 times higher than nominal; the net result is that the faint emission from the galaxy is 
undetectable. The NICMOS frames of NGC1291 missed the galaxy because of guide star 



^The extinction curve k(A) is defined through the foUowing equation: 7 obs{X)^^ int{X) IQ-^-'^^We(b-v) ^ 
where Yobs and Yint are the observed and intrinsic fluxes, and E(B— V) is the color excess. 
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problems. The F187N images of NGC4631 show a faint flat-field imprint (generally a sign 
of residual pedestal) that, coupled with the large dynamical range of the emission from this 
edge-on galaxy, produces a very uneven background. For NGC3034 (M82), problems related 
to non-linearity corrections and saturation for this bright target exist for the NICMOS, IRAC, 
and MIPS images, making photometry in the center of this object highly unreliable at the 
present time. 

The HST archive was also mined for Ha images for those cases where (a) coverage 
was similar between NIC3 and optical images, and (b) the narrowband filter provides a 
better rejection of the [Nil] emission line than the ground-based images. WFPC2 images 
that met these criteria were available for NGC1512, NGC4736, NGC4826, and NGC5055. 
The line emission was observed through the narrowband filters F656N or F658N, and the 
underlying continuum through F547M, F555W, and/or F814W (equivalent to medium-V, 
V, and I, respectively). For NGC4736, NGC4826, and NGC5055, the [Nil] /Ha values hsted 
in Table [2] come from the comparison of the fluxes in the HST and ground-based (see below) 
narrowband filters; the [Nil] contamination in the HST filters is minimal, and has been used 
to guide our extrapolation of the best nitrogen-to-Ha ratio to attribute to each galaxy. This 
value has been used for those areas in the ground-based images not covered by the HST. 



3.3. Ground— based Optical Imaging Data 



R-band and Ha-centered narrowband images were obtained for most of the galaxies 
as part of the SINGS ancill ary data prograrn, eit her at the 2.1-m KPNO telescope or at 
the 1.5-m CTIO telescope (IKennicutt et al.ll2003l ). Exposure times were typically around 
1800 s for the narrowband filters, and a few hundred seconds for R. Standard reduction 
procedures were applied to all the images. Standard stars observations were obtained during 
each observing run to derive photometric calibrations. 

The rescaled broadband images were subtracted from the narrowband images to obtain 
emission-line-only images. The [Nil] contamination within the filter bandpass is removed us- 
i ng [Nil] /Ha values me asured either from large-aperture (~50") SINGS optical spectroscopy 
( iMoustakas et al.ll2007l ) or retrieved from the hterature (Table |2]), and accounting for changes 
in the filter transmission between the wavelengths of Ha and the two [Nil] emission lines. 
High metallicity galaxies for which [Nil] /Ha ratios are not available from either source, or 
cases w hich have optical spect ra dominated by a central non-thermal source (Seyfert 2 or 
LINER. IMoustakas et al.l 120071 ) are assumed to have [NII]/Ha ~0.5. Within each galaxy, a 
constant [Nil] /Ha is adopted, although th e ratio can change signific antly from individual 
HII regions to the more diffuse component (iHoopes &: Walterbosll2003l ). Radial variations of 
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[Nil] /Ha within a galaxy are less of a concern here, as only the central region of each galaxy 
is imaged. 

Typical 1 a sensitivity limits of the final Ha images are 1-2x10^^^ erg s~^ cm~^ arcsec^^, 
i.e., they are a factor 3-10 deeper than the Paa images. This, coupled with the fact that 
the Ha is, intrinsically, about 8 times brighter than Paa, implies that our Ha measurements 
will have higher signal-to-noise ratio than the Paa ones for Ay <4 mag. 

Narrowband and R-band images of DDO053, MSlDwB, Holmberg9, and NGC4625 were 
obtained using a CC D imager on the Stewar d Observatory Bok 2.3 m telescope, as part of 
the IIHUGS project (IKennicutt et al.ll2007bl ). Narrowband and R-band images of NGC5408 
were obtained at the CTIO 0.9-m telescope, also as part of the IIHUGS project. Images 
were taken using a 70 A narrowband filter centered at 6580 A and an R-band filter and a 
Loral 2kx2k CCD detector. Exposure times were 1000 s in Ha and 200 s in R, and reach 
comparable depth to the KPNO images because of the high throughputs of the filter and 
the CCD detector. Data reduction followed similar procedures as described above. 

Ground-based Ha images for NGC3627, NGC4736, NGC4826 and NGC5055 were pro- 
vided by the SONG collaboration Jsheth et al.ll2002l : iHelfer et al.ll2003h . as SINGS did not 
repeat these observations. The data were obtained at the KPNO 0.9-m telescope, with an 
observing strategy and filter selection similar to those of SINGS. The main difference between 
the SINGS and SONG Ha images is the total exposure time (and the depth of the images), 
being in the latter case 3-5 times shorter than in the former. For this reason, the ground- 
based SONG images were used in conjunction with the HST Ha images for photometric 
measurements in NGC4736, NGC4826, and NGC5055. 



4. Photometric Measurements 

4.1. Aperture Photometry 

For each galaxy, the Ha, stellar-continuum-subtracted 8 /xm, and 24 /xm images were 
registered to the same coordinate system of the Paa image, before performing measurements. 
Photometric measurements at all four wavelengths of local 24 /xm and Ha peaks were per- 
formed on the common field of view of the four images. Emission peaks at 24 /xm (and 8 /im) 
have generally corresponding Ha peaks; the opposite, however, is not always true, and there 
are some cases of Ha emission peaks without corresponding mid-IR emission. Thus, both 
24 /im and Ha images were used independently to locate local peaks of star formation. 

The size of the aperture used for photometric measurements is dictated by the lowest 
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angular resolution image, the MIPS 24 fim image, with a PSF FWHM~6". We chose 
apertures with 13" diameter as a compromise between the desire to sample the smallest 
possible scale compatible with HII regions and the necessity to have reasonable aperture 
corrections on the photometry (Figure [T]). For the chosen aperture size, corrections to 
infinite aperture are 1.045, 1.05, and 1.67 at 3.6 /im, 8 fim, and 24 /x m, respectively, for 



point sources (SSC IRAC Handbook and MIPS Handbook, respectively; iReach et al. 



Engelbracht et aL 



2007 



Jarrett 



2005 



and Ha images (jCalzetti et al. 2005). 



20061) . and are assumed to be small or negligible in the Pao; 



In the case of the IRAC 3.6 /im and 8 /xm emission, exten ded emission has a different 
aperture correction than point sources. Best current estimates (jjarrettll2006l ) indicate that 
our aperture choice requires an additional correction factor 1.02 at 3.6 /xm and 0.90 at 8 fim, 
for extended sources. As our sources are neither totally extended nor point-like, actual 
aperture corrections are likely to be closer to a value of unity than those reported here. 

The fixed aperture corresponds to different spatial scales in different galaxies, as dis- 
tances between ~0.5 Mpc (spatial scale ~30 pc) and ~20 Mpc (~1.26 kpc) are covered. In 
order to allow comparison among luminosities measured over areas that differ by a factor 
as much as ~40 (for the typical distance range 3-20 Mpc), we report all measurements as 
luminosities per unit of physical area (luminosity surface density, LSD) Spaa, ^Ha, Ss/^m, dust, 
and S24^m, in units of erg s~^ kpc~^. Luminosities at mid-infrared wavelengths are expressed 
as z/L(z/). 

The use of luminosity surface densities removes most dependence of our measurements 
with distance, as the LSDs are, for our purposes, equivalent to fluxes. Notable exceptions 
are the cases where the area covered by our aperture contains only one HII region, with 
intrinsic size smaller than our adopted fixed aperture's size; in these cases the LSDs will 
be artificially decreased by the larger area of the aperture relative to the values they would 
have if we selected apertures matched to the intrinsic size of each HII region/complex. The 
latter choice is not easily applicable to our sample due to the angular resolution limitations 
of some of the data. Furthermore, we will see in section 6 that this effect does not appear 
to have an important impact on our results. 

Photometry for a total of 220 separate HII knots is obtained in the 33 galaxies. Of 
these, 179 are in the 23 high metallicity galaxies, including 11 non-thermal nuclei (Seyfert 2 
or LINERs as retrieved from NED|f|; no aperture was laid on top of the active nucleus of 



•^The exact classification of galactic nuclei is beyond the scope of the present work; we restrict oursclf to 
well-known non-thermal sources as described in the literature, as these are the sources that most deviate 
from the general trends described in the following sections. 
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the edge-on galaxy NGC5866). In the five medium metallicity and five low metallicity 
galaxies, 22 and 19 regions are measured, respectively, including 4 regions (one each in 
IC2574, Holmberg IX, MSlDwB, and NGC6822) that are strongly emitting in the mid- 
infrared, but are undetected in both our Paa and Ha data. These line-undetected objects are 
detected in the optical continuum bands and are extended; t hus they are lik e ly ba ckground 



sources. Heavily obscured sources, like those discussed in iPrescott et al.l (120071 ) . should 
represent about 3% of the 24 //m sources, but we find none; we attribute this lack of heavily 
obscured sources in our sample to the small spatial region subtended by the NICMOS FOV 
within each galaxy. The 11 non-thermal sources and the 4 background sources (Figures [3H1]) 
will be excluded from all subsequent statistical analysis. 

Crowding of emission peaks within each frame prevents the use of 'annuli' around in- 
dividual apertures to perform background subtraction from the photometric measurements. 
B ackground rern o val is thus achieved by subtracting a mode from each frame, as described 



m 



Calzetti et all (j2005|). 



'Integrated' values of Ha, Paa, 8 /xm, and 24 /im luminosity surface density are also 
derived for each galaxy within the area imaged by the NICM0S/NIC3 camera. These 
integrated values are therefore the LSD of each galaxy within the central ~50", except 
for NGC5194, where the central ~144" are measured (Table [2]). The integrated values mix 
the emission from the star forming regions (measured with the smaller apertures) with areas 
of little or no star formation, thus providing some insights into the impact of the complex 
galactic environment on SFR calibrations. 



4.2. Uncertainties in the Photometric Measurements 

The uncertainties assigned to the photometric values at each wavelength and for each 
galaxy are the quadrature combination of four contributions: Poisson noise, variance of 
the background, photometric calibration uncertainties, and variations from potential mis- 
registration of the multiwavelength images. The variance on the image background is de- 
rived in each case from the original-pixel-size images. The impact of potential background 
under- or over-subtractions varies from galaxy to galaxy, and also depends on the relative 
brightness of the background and the sour ces. The effect of pot ential misregistrations have 



been evaluated for the case of NGC5194 by lCalzetti et al.l (120051 ). Because of the large aper- 



tures employed for our photometry, this contribution is either small (a few % of the total 
uncertainty) or negligible. 

For the Spitzer 8 fim and 24 fj,m images, calibration uncertainties are around 3% and 4%, 
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respectively (IReach et al.ll2005l : lEngelbracht et al.l 120071 ). This, added in quadrature to the 



other uncertainties, produces overall uncertainties in the measurements that range between 
15% and a factor of two, with the median value being around 22%. The superposition of the 
PSF wings in adjacent apertures produces an additional effect in the 24 nm measurements , 



that is evaluated and removed on a case-by-case basis (see example in lCalzetti et al.ll2005l ). 



For the HST images, photometric calibrations are generally accurate to within ~5%, 
for narrowband filters. The faintness of the Paa emission, and therefore the impact of 
the background variance and stellar continuum subtraction is what mostly dominates the 
photometric uncertainty on the Paa emission line measurements, with values between 15% 
and a factor of roughly 2, with a median value of 60%. For the extinction-corrected Paa 
luminosities, the uncertainty on the attenuation Ay increases the Paa uncertainty by a factor 
of 1.22. 

For the ground-based Ha images, which are the deepest images in our set, the main 
sources of uncertainty are: photometric calibrations, stellar continuum subtraction, and the 
correction for the [Nil] contribution to the flux in the narrow-band filter. These translate 
into uncertainties in the final photometric values between 10% and 50% (with occasional 
factor-of-2 uncertainty). The median uncertainty for the Ha luminosities is 20%. Although 
less deep, the HST Ha images are characterized by more stable photometry, better continuum 
subtraction, and smaller [NH] contamination; uncertainties on the final luminosities are in 
the range 5%-10%. 

For a few of the galaxies of Table 2, some special circumstances are present or special 
treatment was required. For NGC2841, the very faint line emission produces large, and 
highly uncertain. Ay values. For NGC5033, no Ha image is available; the uncorrected Paa 
can be up to 70% underestimated for the largest Ay measured in our sample (Ay ~4 mag), 
and, therefore, this galaxy is excluded from all fits reported below. 

In Holmberg IX, Ha emission is detected in two of the three selected regions; for one 
of these two regions, 24 fim emission is also detected, at the ~2.5 a level. A strong 24 fim 
detection is present in the third region, together with the only 8 fim detection in the field; 
because of the absence of hydrogen line emission and of the extended nature of the broad 
band emission, this source is identified with one of the background sources discussed in 
section 4.1. For the two regions with Ha emission, only upper limits can be derived for the 
Paa and 8 /im emission. The presence of Ha emission provides a lower limit to the Paa line 
intensity for the zero extinction case (after including the uncertainty on the Ha measurement 
itself). We have taken the range between this lower limit and the upper limit measured from 
the HST/NICMOS images to be our fiducial range of values for Paa, and therefore we report 
the middle values (in logarithmic scale) as measurements, rather than use the actual upper 



- 16 - 



limits. 

In NGC5408, the brightest, and most extended, hne-emitting region is only partially 
imaged by NICMOS. The Paa image is therefore used only to derive a typical Ay value for 
the region, using small-aperture photometry and the matching Ha measurements. The A^^ 
value derived in this way is then applied to the Ha emission of the entire, extended, region, 
for which a larger-than-nominal, 17.1" diameter, aperture is used, not only for Ha, but also 
for the 8 /xm and 24 /xm emission. The other two regions in this galaxy are treated with the 
nominal procedure described in section 4.1. 



5. Starburst Galcixies 



Our baseline sample of 220 HII knots is augmented with 10 local low-met allicity star - 
burst galaxies and 24 LIRGs from lEngelbracht et al.l (120051 ) and lAlonso-Herrero et al.l (120061 ). 
respectively, in order to verify that trends and correlations observed for star-forming regions 
within galaxies can also be applied to galactic-scale (~kpc) star formation. In this context, 
starbursts are defined as galaxies with a central, connected star forming region whose energy 
dominates the light output in the wavebands of interest. 

The low-metallicity starbursts and the LIRGs also expand the mid-IR and line emission 
LSD parameter space of the low- and high-metallicity HII knots, respectively, by more than 
an order of magnitude at the high end. 



5.1. Low— Metallicity Starburst Galaxies 

As part of the HST/NICMOS SNAP-9360, about 40 nearby starburst galaxies were 
obser ved. Of these, 13 also hav e Spitzer imaging as part of the MIPS and IRS GTO observa- 



tions (lEngelbracht et al.ll2005l ). The main characteristics and measurements for 10 of these 
galaxies are listed in Tables [3] and HI The three remaining galaxies, NGC3079, NGC3628, and 
NGC4861, are omitted from the present analysis for the following reasons. For NGC4861, 
the HST/NICMOS pointing targeted the relatively quiescent center of this galaxy, rather 
than the peripheral giant HII region. The other two galaxies, NGC3079 and NGC3628, have 
extended optical line and mid-IR emission: about 40% and 60% of the emission is outside 
of the field-of-view imaged by HST/NICMOS; corrections for the fraction of light in the 
Paa line outside of the observed frame would be thus substantially larger than the typical 
uncertainties in the measurements. 



The data for the galaxies in Table [3] were reduced in the same fashion as the SINGS 
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galaxies discussed in sections 2-3. In particular, the HST/NICMOS images, which are 
presented here for the first time, were treated following the same procedure as section 3.2. 
The main difference between the HII knots in the SINGS galaxies and the local starbursts 
is in the photometry: integrated flux values encompassing the entire central starhurst (the 
dominant source of emission at the wavelengths of interest) a re derived for the la t ter sa mple. 
The integrated measurements at 8 /xm and 24 /xm are from lEngelbracht et al.l (120051 ). and 
are reported in Table HI 

The Paa measurements (Table HI) are performed using the aperture sizes listed in Table El 
and are corrected for the Galactic foreground extinction (fourth column of Table [3]), but not 
for internal extinction. We expect the internal extinction to represent a small effect on the 
Paa flux in these mostly low metallicity gal axies (compar e with Figure [2]). An exception may 
be represented by SBS0335-052, for which iHouck et all (120041 ) measure Ag 7 ~0.5 mag. 
If the region of silicate absorption is coincident with the region of line emission, this would 
correspond to P^paa ~2 m ag. Given th e unce rtainty in the spatial co-location of the dust- 
hidden source detected by iHouck et al.l (120041 ) and the main source(s) of the Paa emission 
and the fact that the introduction of an extinction correction for one of the galaxies does 
not impact our conclusions, we do not perform the correction. 



5.2. Luminous Infrared Galaxies 



HST/NICMOS Paa data and extinction corrections, as well as information on the physi- 
ca l extent of the star forming ar ea for each of the 24 LIRGs used in this analysis, are presented 



m 



Alonso-Herrero et al.l (120061 ) ; the reader is referred to that work for details. Infrared mea- 



sure ments at 25 jUm from IRAS and distances for each galaxy are from ISanders et al.l (120031 ) 
and ISurace. Sanders fc Mazzarellal (120041 ). At the time of this writing, no 8 /xm emission 
measurements are available for these ga laxies. The LIRGs' metallici ties are characteristic of 
our high-metallicity HII knots sample (lAlonso-Herrero et al.ll2006l ). Photometry for these 
galaxies, as in the case of the local starbursts (section 5.1), includes the entire line-emitting 
and IR-emitting galactic region, thus the measurements are integrated galaxy values. 



6. Analysis and Results 

Photometric measurements for the 220 HII knots, the local low-metallicity starbursts, 
and the LIRGs are shown in Figures [SHU where the infrared LSD in the two mid-IR wave- 
bands is shown as a function of the extinction-corrected Paa LSD. 
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One characteristic immediately apparent in Figures [3H1] is the overall correlation between 
the infrared LSDs and the Paa LSD (panels [a]), although the scatter is non negligible in 
both cases (panels [b]). The correlations appear especially significant for the high metallicity 
HII knots (the most numerous subsample among those under analysis here), and span a 
little over two orders of magnitude in Paa LSD. Bi-linear least-square fits through the high- 
metalhcity data points yield: 



Log Ss ^m, dust = (0.94 ± 0.02) Log Spaa,corr + (4.80 ± 0.85), 



(2) 



Log S24 = (1.23 ± 0.03) Log Spaa,corr - (6.88 ± 0.97), (3) 

where Spaa,corr is the extinction-corrected Paa LSD. Equation 3 accounts effectively for the 
trend of the LIRGs, although these data were not used in the fitting procedure. 

The scatter of the datapoints about the best fit lines of equations 2-3 are approximately 
the same, with cr=0.3 dex (panels (b) of Figure [3H1]) . Thus the 1 cr scatter is about a factor 
of 2 for the high metallicity regions. 

A potential source of bias in equations 2 and 3 is the large range of distances that our 
sample covers, about a factor of 6 for the high metallicity galaxies. Our fixed photometric 
aperture of 13" diameter thus probes regions that are about 30 times different in area between 
the nearest and the farthest targets in the high metallicity subsample, i.e., from 0.04 kpc^ 
at 3.5 Mpc to 1.12 kpc^ at 17 Mpc (for the most distant galaxy in our sample, NGC4125, 
located at 21 Mpc, only the central Sy2 nucleus is detected and is excluded from the analysis). 
Although we remove the background from each photometric measurement, uncertainties in 
this subtraction will affect the farthest targets more strongly than the closest ones, if HII 
regions/complexes have constant sizes of ^5^100-200 pc. Furthermore, we may expect that 
our fixed aperture photometry may dilute the LSDs of the more distant regions, for the 
extreme hypothesis that only one HII region is contained in each aperture. We have tested 
the impact of these effects by looking at the distribution of the ratios Ss /^m, dust/Spaa,corr 
and S24 fj.m/Spaa,corr SiS & fuuctiou of galaxy distance (Figure [5]). For the high metallicity 
subsample, non-parametric (both Spearman and Kendall) tests show that the data are 
uncorrected with the galaxy's distance, suggesting that there is no obvious bias in our 
analysis. 

Both the 8 fim and 24 /zm LSDs of medium and low metallicity regions are deficient 
relative to the extrapolation of the best fit lines for the high metallicity regions (Figure [3] 
andHj). The deficiency is far more pronounce d in the case of Ss/,,m, rf,;,.st i a fact already noted 



in a number of previous investigations fe.g.. lEngelbracht et al. 2005; Galliano et al. 2005 



Hogg et al.ll2005 



Rosenberg et al.ll2006l : iDraine et al.ll2007l ). A potential source of concern in 



this case is that the high metallicity subsample has a higher mean distance than the medium 
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and low metallicity ones (Figure [5]). iHelou et al.l (12004 ) have shown that the 8 fim emission 
is brighter at the edges of an HII region (i.e., in the PDR) than at its center. Our fixed 
aperture photometry could therefore underestimate the 8 fxm flux from the low metallicity 
regions, if the apertures are not large enough to sample the entire area surrounding the 
HII knot. However, Figure [5] shows that the 8 /zm emission is deficient in the medium 
and low metallicity subsamples relative to the high metallicity one even when galaxies at 
comparable distances are considered. The only potential exception is NGC 6822, the closest 
galaxy to the Milky Way in our sample, which, at a distance of 0.47 Mpc, could suffer from 
the effect of having too a small aperture applied to the 8 /im emission measurements; indeed 
its mean value is lower (although not statistically significantly) than the average of the other 
data in the same metallicity bin. 

The trend of the Ss^m, dust/S24:fim ratio as a function of Spaa,corj- (Figure ED highlights 
the decrease of the 8 /im LSD for decreasing metallicity, and also shows that the effect is 
independent of the number of ionizing photons in the region. The latter suggests that: (1) 
our aperture sizes are large enough to encompass both the HII regions and the surrounding 
PDRs, as noted above; and (2) in these large regions the dependence of the 8 /xm-to-24 //m 
ratio on the luminosity surface density of the HII region/complex that heats the dust is a 
small effect relative to the effect of metallicity. The decrease of the 8 /xm to 24 /im LSD ratio 
as a function of increasing Paa LSD in the high metallicity points (i.e., at roughly constant 
metallicity) indicates that the component of thermal equilibrium dust contributing to the 



Helou 


1986: 


Draine & Li 


2006) 



contribution may also come from an increased destruc tion rate of the 8 fim dust emission 



for increasing starlight intensity (IBoulanger et al.lll988l ). 



As suggested by lKennicutt et al.l (l2007al ). the combination of measurements at Ha and 
24 yum can provide insights into both the unobscured and obscured regions of star formation. 
We have combined linearly the observed Ha and 24 fim LSDs and scaled them to the Paa 
LSD. The best fit line through the data is: 



Log 



Paa,corr 



(0.98 ± 0.02) Log {aSHa,obs + hS^A ^m) + (0.73 ± 0.93) 



(4) 



where a is the intrinsic Paa/Ha ratio, thus is dictated by atomic physics and is only mod- 
erately dependent on metallicity (a=0.128, 0.118, and 0.114 for the high, medium, and low 
metalhcity data, respectively; see section 3.2). The coefficient h for the 24 /xm LSD has been 
empirically rescaled to bring the sum of the optical and IR LSDs in agreement with the Paa 
one (6=0.0040, 0.0037, and 0.0036 for the high, medium, and low metallicity datapoints, 
respectively; Figure [7]). The best fit from equation 4 gives 6/a=0.031±0.006, and this ratio 
is independent of metallicity. Equation 4 is, within the uncertainties, consistent with a linear 
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relation with null intercept between the two quantities, as expected if the right-hand-side 
expression is a measure of the ionizing photon rate, like Spaa,corr- The linearity of the rela- 
tion is by construction, as the requirement is to approach unity as much as possible for all 
the combined data, but the null intercept has not been fixed a priori; furthermore, the ratio 
b/a was left as a free parameter in the analysis, and its constant value is a result (not an 
input). 

Interestingly, the high metallicity datapoints show approximately the same dispersion 
around the mean trend of equation 4 as they do for equations 2 and 3, with a 1 a ~0.3 dex. 
In the case of the combined optical/mid-IR, the dispersion is the same whether the high 
metallicity datapoints alone or all datapoints are included in the statistical analysis (panel (b) 
of Figure [7]). Conversely, for the two mid-IR LSDs the dispersion is measured for the high 
metallicity datapoints only, and increases substantially (on one side) when the medium and 
low metallicity datapoints are included in the statistics (panels (b) of Figures [3] and H]). 
These considerations do not include the LIRGs, that in Figure [7] show evidence of having 
higher combined optical/mid-IR LSDs than inferred from the extrapolation of equation 4. 
A possible explanation for this effect will be discussed in Section 7. 



As already discussed in iKennicutt et al.l (l2007al ). the sum on the right-hand-side of 



equation 4 can be interpreted as a representation of the dust extinction corrected Ha lumi- 
nosity or LSD. as: 

SHa,corr = SHa,obs + (0.031 ± 0.006)S'24 ^m- (5) 

The proportionality c oefficient for the 24 /iin luminosity is ~20% smaller than that derived 



for NGC5194 alone (IKennicutt et al.l l2007al ). which is within the 1 a uncertainty. This 



small difference is likely due to the larger variety of galaxies used in the present work which 
provides a dynamical range in luminosity surface density about an order of magnitude larger 
than in the NGC5194 case. 

The proportionality coefficient for the 24 fim emission in equations 4 and 5, b/a=0.031, 
is independent of metallicity. This suggests that in the S24/im versus Spaa,corr plane the 
observed deviations of the medium and low metallicity data from the best fit for the high- 
metallicity datapoints are simply due to the progressively lower dust content of the ISM 
for decreasing metallicity (section 7). No other effect beyond the simple increase in the 
medium's transparency is required. Indeed, most of the contribution to SHa,corr comes 
from the observed Ha emission at low Spaa,corr LSDs (low dust systems) and, vice- versa, 
it is mainly contributed by the 24 fim emission at the high LSD end of our sample (dusty 
systems) . 



Similar correlations as those seen for the HII knots within galaxies exist between the 
integrated LSDs of the galaxies' centers (section 4.1). Figure [8] shows the trends for the 33 
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star-forming galaxies in our main sample. For the combined optical/mid-IR LSDs, a linear 
fit through the integrated datapoints of the high metallicity galaxies are consistent, within 
1 0", with the best fit lines through the individual HII knots, both in slope and intercept 
(third panel of Figure [H]). For the 24 /im LSD, the slope of the linear fit is consistent (again 
within 1 a) with that of the individual HII knots, and the intercept is consistent (within 
0.1 a) with the value expected by simply rescaling the HII knots' mean LSD for the larger 
area used in the integrated measures. The results for both the 24 /xm and the combined 
optical/mid-IR integrated measures suggest that within the central areas covered by the 
NICMOS observations any diffuse 24 fim emission contributing to the measured LSD is 
matched by diffuse Paa LSD with comparable intensity. This, of course, does not mean that 
diffuse 24 /im emissi on is not present ; indeed, such diffuse emission has been observed in 



the SINGS galaxies (iDale et al.ll2006l ). Our result simply implies that such diffuse 24 /xm 
emission traces the diffuse ionized emission, at least within the central galaxy regions sampled 
by our data. 

A more complicated scenario appears for the 8 /im LSD: a best fit line through the 
high metallicity integrated regions produces a higher slope (1.16±0.09) than derived for the 
individual HII knots. The difference is marginally significant (2.2 a), but implies that the 
8 /im LSD is higher by about a factor of 2 over what is expected from a simple rescaling of 
areas at the high luminosity encj^. A visual inspection of the images shows that the galaxies 
with low 8 /im LSDs generally have line and mid-IR emission which is centrally concentrated 
or coming from thin, almost edge-on, disks or annuli located in the central 50", while at the 
high 8 /im LSD end galaxies tend to have a more homogeneous distribution of HII knots. 



7. Comparison with Models 

To help clarify the nature of some of the characteristics of the observed correlations, 
this section is devoted to the comparison of our data with simple models that exploit the 
energy balance between the stellar light absorbed by dust at UV, optical and near-infrared 
wavelengths, and the light emitted by the dust in the mid- and far-infrared. The details of 
the models are presented in the Appendix, in addition to a discussion on limitations to their 
use and applicability. Here we provide a brief summary of those models. 

The basic approach adopts a range of plausible stellar populations for our HII knots 
(and starburst galaxies), in terms of star formation histories, stellar population ages, and 



''For the 'integrated' diffuse emission, tlie extended source aperture correction provided by ljarretti (200] 
has been used. 
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metallicities (2005 update of StarburstQ^lfl, Leitherer et al. 1999 ). Simple assumpti 



tions are 



then dust-attenuated accordinEr to empirical recipes ( 


Ca 


zetti. Kinnev & Storchi-Beremann 


1994; 


Meurer. Heckman & Calzetti 


1999: 


Calzetti et al. 


2000; 


Calzetti 


2001 


) to provide a 



'predicted' infrared emission, S/h- As the stellar populations probed in our analysis range 
from groupings of a few to several HII regions for the HII knots to populations with extended 
star formation histories in the case of starbursts and LIRGs, both instantaneous bursts and 
constant star formation populations are included. The total infrared emission will, in gen- 
eral, depend not only on the adopted stellar population, but also on the extinction curve and 
the dust geometry. Sin ce for th e last t wo parameters, we make a simplifying assumption and 
use the prescription of ICalzettil (l200ll ): the impact of varying the dust geometry is discussed 
in section A. 4. For the spectral e nergy distribution (SED) of the infrared emission, S/i^, we 
adopt the model of iDraine &: Lil (120061 ). according to which the fraction of IR power emerg- 
ing in the IRAC 8 /zm and MIPS 24 /im bands is a function of the starlight intensity. We 
determine (section A. 2) the range of starlight intensities corresponding to the model stellar 
populations we are considering, so to obtain a direct correlation between the Paa LSD and 
the fraction of IR light emerging in the two mid-IR bands. Since our HII knots follow the well 
known correlation between SFR and ext i nction (Section A.l and Wang &: Heckma"nl 19961: 
Heckman et al.lll998l : iHopkins et al.ll200ll : ICalzettill200ll : iMoustakas. Kennicutt &: Tremonti 
20061 ). which we parametrize as a relation between color excess E(B— V) (section 3.2) and 
the ionizing photon rate per unit area Ejo„, we use this relation to link the stellar population 
models to the dust attenuation model, and eliminate one degree of freedom in our models. 
Model parameters that we allow to vary are the star formation history of the stellar popu- 
lations (bursts or constant star formation), their age (0-10 Myr for i nstantanous bur s ts, th e 
range chosen to ensure presence of significant ionizing photon rate, iLeitherer et al.l (119991 ): 
6-100 Myr for constant star formation), the mass (10^-10*^ Mq) or SFR (4x10^^-4 Mq yr^^) 
of the stellar cluster(s) associated with the HII knot or starburst galaxy, and the metallicity 
of both the population and the interstellar medium (0.1-1 zj^. Figures [914TT] show the basic 
results from the comparison between the models described so far and our data for the 8 /xm, 
24 /im, and combined optical/mi d-IR emission from HII knots and star-forming galaxies. 

The larger-than-unity slope of the 24 yum versus Paa LSD (in log-log scale. Fig- 
ure [9]) is a natural outcome of the models in the high luminosity surface density regime, 
Log(Spaa,ofes)>39, and is an effect of the 'hotter' IR SEDs for increasing starlight intensity. 



^http://www. stsci.edu/science/starburst99/ 

^We adopt the oxygen abundance 12+log(0/H)=8.7 as solar metallicity value (jAllende Prieto et al ]l200ll) . 
which we take here as representative of our high-metallicity HII knots. 
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In other words, regions with higher Paa LSD emit proportionally more of their infrared 
energy into the 24 fim band, because the peak of the IR SED m oves towards shorte r wave- 



lengths (higher 'effective' dust temperatures, see Appendix and iDraine &: Lil (120061 ) ). The 
models also predict a slightly larger than unity value for the slope of the 8 fim LSD correla- 
tion, which is steeper than that of the HIT knot data (Figure [T0|) . but is roughly consistent 
with the slope of the integrated measures. 

The models account well for the linear relation of the combined optical/mid-IR LSD 
with the Paa LSD (Figure [TTi). for luminosity surface densities Spaa,corr <10^° erg s~^ kpc~^. 
At high luminosity surface density, the models for the combined LSDs depart from a linear 
relationship, as increased starlight intensities are expected to raise the temperatures of the 
larger grains so that the fraction of the absorbed energy re-radiated at 24 /^m (which is, at 
these high LSDs, the dominant contribution to equation 5) increases. The LIRGs data, that 
populate the high LSD regime in our plot, do indeed confirm observationally the deviation 
from the extrapolation of the best fit line; they show a steeper-than-one slope, in qualitative 
agreement with the models' expectations (Figure [TTIl . 

At the high luminosity end (LIRGs and brighter), an additional effect that can con- 
tribute to the deviation from the slope of unity observed in Figure [11] and the steeper-than — 
unity slope of Figure [9] is the competition between the dust and the gas for the absorption 
of some of the ionizing photons. In the high luminosity regime, st ar formation occurs i n 



environments of increasing density, e.g., ultracompact HII regions (iRigby &: Riekd l2004j ) . 
and the dust absorbs the ionizing photons before they can excite the gas. In this regime, 
standard extinction-correction methods become progressively less effective at recovering the 
intrinsic Pact emission, and will produce an underestimate of the hydrogen emission line LSD 
at constant 24 fim LSD (section A4). The impact of this effect on our data is unclear (and 
currently not included in our models), although it may be relatively small as the bulk of the 
observed trends is fully accounted for by our baseline model. 

Instantaneous burst populations and constant star formation populations produce mostly 
degenerate models for all three mid-IR quantities (Figures [9]- [TTl) . A young, 4 Myr old, in- 
stantaneous burst population in the mass range ~10^-10'^ Mq provides similar model lines 
as a constant star formation model forming stars since 100 Myr and with SFR in the range 
4xlO-M M© yr-^ 

However, even the high-metallicity HII knots in Figures [914T0] show a fairly large disper- 
sion around the mean trends described above, with a clear increase of the dispersion around 
the mean Sg ^m, dust and S24 values for Spaa,corr <10^^ erg s^^ kpc~^. Furthermore, in this 
Paa LSD regime, most of the 8 /xm and 24 fim emission from the high-metallicity HII knots 
is located above the baseline model lines, i.e., the models underpredict the mean values of 
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the mid-IR emission (Figures [9]- [T0l) . The 'downward' curvature of the models is a direct 
product of the increasing transparency of the interstellar medium for decreasing ionizing 
photon rate density and, from equation A2, decreasing dust amount. With a more trans- 
parent medium, proportionally less IR radiation is produced. The medium is still thick to 
Lyman continuum photons, and the ionized hydrogen emission lines are still a measure of 
the total number of ionizing photons in the region. An additional parameter is thus required 
to account for both the large scatter of the datapoints around the mean trends and the large 
number of high— metallicity datapoints above the model lines for the Sg^^. dust and 824/^™ 

LSD 

plots. This second parameter appears to be the age of the stellar population. Ageing bursts 
between ~0.01 Myr and ~8 Myr produce a decreasing number of ionizing photons, while 
at the same time remaining luminous at UV-optical wavelengths (the major contributors 
to the IR emission). Figures lUHTT] show that the 'flaring' of the high-metallicity HII knots 
datapoints around the mean value for decreasing Paa LSD is compatible with the 'flaring' 
of the ageing burst models. Such ageing populations can also account for the data points 
above the mean trends in Figures [9] and [TOl 

The presence of ageing bursts is a sufficient (and physically expected), but not a neces- 
sary, condition to account for the dispersion in the data. As briefly discussed in the Appendix 
(section A. 4), different assumptions from our default one about the average dust geometry 
can also produce a higher mid-IR emission than our flducial model lines. For instance, pres- 
ence of ultracompact HII regions within our HII knots will produce higher IR emission at flxed 
Spaa,corr than cxpcctcd from the models. This is a consequence of the higher opacity of such 
regions, for which the use of the Ha/Paa ratio to recover the intrinsic line fluxes will l ead to 



an underestimate of the intrinsic Paa luminosity in the region. Recently, iDale et al.l (120061 ) 
have shown that for local star-forming galaxies the UV/IR ratio is heavily determined by the 
morphology of the 24 /xm dust emission, in particular by the 'dumpiness' of such emission, 
which therefore determines the escape fraction of UV photons from star-forming regions. 



extinction and attenuation used in 


;he present work ( 


Calzetti, Kinnev & Storchi-Beremann 


1994; 


Meurer. Heckman & Calzetti 


1999: 


Calzetti 


2001 


)■ 



For the combined optical/mid-IR LSD, the models are degenerate as a function of 
metallicity (Figure [TT]) . This is not surprising if the main driver of the discrepancy between 
the high and low metallicity S24^m at flxed Paa LSD is the larger medium transparency, 
i.e., lower dust column density, in the lower metallicity data (equations A2 and A4). This 
is indeed the case (Figure [H]): the separation at low Paa LSD between the solar metallicity 
and the 1/lOth solar metallicity model lines is mostly due to the metallicity scaling factor in 
equations A2 and A4, and, to a much smaller extent, to the difference in metallicity of the 
two stellar populations. The 1/lOth metallicity model line in Figure [9] provides the lower 
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envelope to the data; most of the galaxies in our sample are above 1/lOth solar in metallicity, 
and thus are expected to lie above this model line. 

This result lends crede nce to the use of a combination of S24^m and S^a.obs (equation 5 



and iKennicutt et al.ll2007al ) as an effective tool for measuring the ionizing photon rates, and, 
ultimately, SFRs, at least up to Paa LSDs f^ilO^^-lO^^ erg s~^ kpc~^. In this framework, 
S24/^m probes the obscured star formation, and the only metallicity effects are those induced 
by reduced opacity; conversely, SHa,obs probes that part of the star formation unabsorbed by 
the dust, independent of the gas metallicity. The behavior of the models in Figure [TT] shows 
little difference between different parameters choices, at least within our data uncertainties, 
and they reproduce the main trend of the data reasonably well. 

The discrepancy observed between the high metallicity and low metallicity Sg^^, dust 
data at fixed Paa LSD re quires one addition al ingredient, together with the increased trans- 



parency of the medium. iDraine fc Lil (120061 ) have suggested that the fraction of low-mass 
PAH molecules present in the dust mixture decreases for decreasing metallicity. In the Ap- 
pendix, we show that the two ingredients (increased medium transparency and decrease of 
low-mass PAH molecule fraction) provide comparable contributions to the depression of the 
8 /im emission, and the two together produce the expected lower envelope to the datapoints 
in Figure [TOl 



8. Discussion 



The scope of this study has been to investigate the extent of the regime of applicability of 
mid-IR emission as a SFR tracer, to use models to reproduce the main characteristics of the 
data, and to investigate reasons for any limitation we have encountered. The general trend of 
mid-infrared luminosity surface densities to correlate with the ionizing photon rates or with 



suits, see 


Roussel et al. 


2001 


; I 


^orster Schreiber et al. 


2004; 


Boselli. Leaueux & Gavazzi 


2004 


Calzetti et al. 


2OO5I; 
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2005: 
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8.1. The Combined Optical/Mid-IR SFR Indicator 

Of the three indicators investigated here, the linear combination of the observed Ha 
and the 24 /xm emission is the one most tightly correlated with the extinction-corrected 
Paa emission. The linear relation between the combined optical/mid-IR emission and the 
SFR as traced by Spaa,corr is common to all galaxies investigated, independent of their 
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metallicity. The most straightforward interpretation (iKennicutt et aLl l2007al ) is that the 
24 fim emission traces the dust-obscured star-formation, while the observed Ha emission 
traces the unobscured one. The combination of the two, thus, recovers all the star formation 
in a region. This interpretation is confirmed by the models investigated in the previous 
section, which also suggest the trend to be relatively independent of the characteristics of 
the underlying star-forming population. The deviations from the linear relation (i.e., from 
a slope of 1 in log-log scale. Figure [TTj) observed at luminosities larger than Spaa,corr >10^°- 
]^q4i gj^g g-i ]jpc-2 gj-Q g^jgQ consistent with the models' expectations: as the 24 /im emission 
starts dominating the luminosity budget, the same physical mechanism producing the S24 fim^ 
versus-Spaa^corr trend also produces that of the combined optical/mid-IR indicator. We 
speculate that this mechanism (see below) is the emission from grains with approximately 
steady temperatures, rather than transiently-heated grains, that come into play at high dust 
temperatures, and which causes the ratio of 24 /xm/IR to deviate from a constant value at 
low starlight intensities to one that increases with the starlight intensity. 

Using our baseline best-fitting model of 100 Myr constant SFR, for solar metallicity 
and the stellar initial mass function (IMF) described in section A. 2, the conversion between 
SFR and Ha luminosity is: 

SFR{Mq yr~^) = 5.3 x 10-*'^L{H a) carrier g s"^). (6) 

Variations of ±20% over the constant in this relation are present for younger ages and metal- 
licities dowr i to ~1 /5t h solar. The ~50% difference between the calibration in Equation 6 



and that of iKennicuttI (jl998al ) is mainly due to differences in the stellar IMF assumptions 
(59%), with a small contribution in the opposite direction coming fror n different assump 



tions on the stellar populations (100 Myr in our case versus infinite age in IKennicuttI (Il998al ). 
which gives a 6% decrease to the discrepancy given by the different IMFs). Using equations 5 
and 6: 

SFR{Mq yr~^) = 5.3 x 10-^^[L{Ha)obs + (0.031 ± 0.006)L(24 /im)], (7) 

where the luminosities are in erg s~^, and L(24 fim) is expressed as z/L(z/). This calibration 
does not change if the luminosities are measured over a substantial area of the galaxy (in our 
case the inner ~0.8 to ~5.1 kpc), rather than in smaller regions hugging the HII complexes 
that produce the ionizing radiation (Figures |7] and |H]). However, the potential non-linearity 
at large LSDs is an important caveat. 



8.2. The 24 /im SFR Indicator 



Conversely, neither the 8 /im emission nor the 24 /xm emission alone are linearly corre- 
lated with the number of ionizing photons that are measured in a region. The non-linearity 
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at high 24 /im luminosity is a dire ct consequence of the increasing dust temperature fo r 
more actively star forming objects (ILi fc Draindl200ll : iDraine &: Lill2006l : iDale et al.ll200ll ): 
higher dust temperatures corres pond to highe r fract ions of the infrared emission emerging at 
mid-IR wavelengths. Following iDraine fc Lil (120061 ). in the regime of low stellar intensities 



(low SFRs in our actively star-forming regions, or roughly Spaa,corr < lO'^^ erg s~^ kpc~^), 
most of the 24 /^m emission comes from single photon transient heating of small grains. In 
this case, the 24 /xm photon flux is directly proportional to the stellar UV photon flux (or 
any other photon capable of single photon heating). Thus, the 24 /im emission counts stel- 
lar UV/optical photons, while the Paa counts the Lyman continuum photons; since there 
is proportionality between the two types of photons, the expectation is for a linear scaling 
between 24 fim emission and Paa emission. However, this regime corresponds to the Paa 
LSD range where the decrease of dust opacity also decreases non-linearly with the amount 
of stellar energy re-processed by dust in the infrared. At high Paa LSDs, hence high stel- 
lar intensities, the dust absorbing most of the stellar photons is warm. There is thus an 
increasing contribution to the 24 /zm emission from larger, warm grains (the Wien side of 
the emission from grains), which leads to a non-linear d ependence of t he 24 /im flux on the 
stellar flux. This is in agreement with the conclusions of ISmith et al.l (120071 ). which observe 
a decrease of the PAH/24 /im luminosity ratio for increasing 24 /im/70 fim luminosity ratio; 
this dependence is highly suggestive of an increasing contribution of warm dust to the 24 fim 
emission. 

The observed non-linearity in the 24 /xm versus Paa relation, S24 ocSp^^ ,,^^^, also 
argues against the case that the high LSDs values measured of our apertures may be due to 
the cumulative contribution of many faint HII regions, rather than a few, increasingly bright 
HII regions. In the case of many faint HII regions (low stellar intensities, and, therefore, 
single photon heating) we should expect the 24 /im LSD to scale linearly with the Paa LSD 
at the high end. The observed non-linear behavior argues in favor of the high luminosity end 
to be contributed mainly by intrinsically bright regions, although presence of apertures with 
many faint HII regions cumulatively giving a high LSD may still be present and contribute 
to the scatter of the datapoints around the mean trend. 

At low metallicities, the deviation from a linear c orrelation is due to lower opacities 
for decreasing metal content, and thus column densiti es (IWalter et al. I l2007h . The effect has 
been well known since the early IRAS observations ( Helou et al. IQSsh : as the metallicity 
decreases, regions become proportionally more transparent and emit less in the infrared as 
a larger fraction of the radiation escapes the area unabsorbed by dust. This accounts for 
the underluminosity of the medium in low metallicity regions in correspondingly intense 
ionizing fields. Using the 24 /xm luminosity as a SFR tracer is thus subject to many caveats, 
including that lower metallicity sources will generally be more transparent than their metal- 
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rich counterparts, and the infrared emission will typically underestimate their SFR by a 
factor ~2-4. 

The non-linear correlation between S24 and '&paa,coTr requires some care for deriving 
SFR calibrations. From equation 3, and using equation 6, we derive a SFR density (SFR 
per unit ) calibration: 



^SFR 



{Mq yr kpc 



-2\ 



1.56 X 10-^^ [^24 {erg s'^ kpc 



-2M0.8104 



In order to derive a calibration for SFRs, we convert our LSDs into luminosities, and the 
resulting best fit through the high-metallicity datapoints produces: 



SFR{Mq yr 



1.27 X 10-^«[L24 Mm (erg s^')] 



-ImO.8850 



(9) 



The exponents of equations 8 and 9 are the same within the 3 a error (the combined 1 a 
uncertainty is 0.03). Equation 9 is closer to a linear relation than equation 8, because we 
add a distance effect when using luminosities (which depend on the distance squared). A 
large sample of regions with comparable distances may be needed to fully sort out intrinsic 
effects from distance-related effects. Both relations are derived from best bi-linear fitting of 
data on HII knots, but their extrapolations account for the observed properties of LIRGs as 
well. Additionally, when considering more extended galactic regions or starburst galaxies, 
equation 8 does not change significantly (Figures H] and [8]). Equations 8 and 9 may thus be 
applicable to galaxies in general whose energy output is dominated by recent star formation. 



Equation 9 is very similar to that of lAlonso-Herrero et al.l (120061 ) . who have derived a 
SFR calibration for the 24 /xm luminosity using a sample of Ultraluminous Infrared Galaxies, 
LIRGs, and NGC5194. The difference in th e calib ration constant between our equation 9 
and the calibration of lAlonso-Herrero et al.l (120061 ) is entirely due to the slight difference 
in exponent betwe en the two relations, and the different SFR-L(Ha) calibrations used here 
and in that work. iPerez-Gonzalez et al.l (120061 ) find a lower exponent, ~0.77, than the one 
in equation 9, about a 4 a difference; however, their result is based on line emitting regions 
in just two galaxies, NGC5194 and NGC3031. 



8.3. The 8 /um Emission 



The analysis of the HII knots in M51 has shown a general, non-linear correlation between 
the 8 fim and the Paa emission (with exponent 0.79. ICalzetti et al.ll2005l ). The present study 
similarly recovers a non-linear behavior for the high-metallicity data, albeit less extreme 
than in the M51 case: Ss^m, dust oc Sp^^^^^^. Our simple 7i=Zq models also predict a non- 
linear correlation between the 8 fim and Paa emission, but with an exponent slightly above 
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unity. Therefore, the gap between observations and expectations is even wider than a simple 
deviation from a hnear correlation; with our uncertainties, the discrepancy is at the 10 a 
level. This level of discrepancy remains unchanged when other uncertainties, e.g., on the dust 
modelling and on the correlation between dust attenuation and number of ionizing photons 
(see Appendix), are included. In contrast, the same simple models are quite successful at 
explaining the observed trend of S24 ^m- We conclude that the 8 /im emission as measured 
within our apertures must include additional contributions that are not included in our 
simplified models. 

Mechanisms that can produce a lower-than-expected slope in a correlation between 
Ss/im, dust and '&paa,corr iucludc the potential contami nation of our measurements by the dif- 



fuse emission from the general galactic radiation field ( Li fc Draine 



2OO2I : Boselli. Lequeux &: Gavazzill2004 ; Peeters. Spoon fc Tielens 



2002 



2004 



Haas. Klaas fc Bianchi 



Wuetal 



20051: 



20051) and/or destruction/fragmen t ation of the 8 /im e r nission carr i ers (jBoulanger et al. 



Mattioda et al. 



1988 



I990I : iHelou. Rvter fc Soifeil Il99ll : iHouck et all I2OO4I : IPetv et all |2005|). In the case of de 



struction or fragmentation of the PAH emitters at 8 /im, the brightest HII regions will show 
a deficiency in the 8 /xm luminosity relative to the fainter regions. In the case where non- 
ionizing populations, as well as ionizing ones, heat the 8 /xm dust carriers, the contribution 
of the former to the Sg^^, dust measurements within our apertures will become proportionally 
larger as the HII regions become fainter (decreasing '&paa,corr), again flattening the observed 
trend. Finally, if the volume filling factor of the 8 /im luminosity originating in the PDRs 
evolves differently from that of the HII regions (or other inhomogen eities in the 8 /im emis- 



sion distribution are pre sent) as the HII region's luminosity increases (IForster Schreiber et al. 



2004j : lHelou et al.ll2004j ). the net result will also be a lower-than-expected exponent between 



dust 



and Sp, 



aa.corr- 



The correlation exponent predicted by models of the 8 /^m LSD versus Paa LSD is 
better matched by the data of large-scale (multiple-kpc) measurements of the centers of 
the high-metallicity galaxies (Figure [8]). There is, however, an offset between models and 
observations, in the sense that the data imply about 65%-100% more 8 /im emission than 
expectations from models. This result argues in favor of one of the mechanisms described 
above: the presence within our apertures of diffuse 8 /im emission unrelated to the current 
star formation. Indeed the result can be explained if stellar populations other than those 
related to the current star formation can heat the carriers of the 8 /im emission, and their 
contribution becomes proportionally larger than that of the ionizing populations as the size 
of the sampled region within each galaxy increases. 



The underluminosity of the 8 /xm emission a t low metallicities (IBoselli. Lequeux fc Gavazzi 



2OO4J : lEngelbracht et al.ll2005l : iHogg et al.ll2005l : iGalhano et al.ll2005l : iRosenberg et al 



2006 
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Madden et al.ll2006l : IWu et al.ll2006l : iDraine et al.l 120071 ) is not a function of the intensity of 
the radiation field (Figure [6]). Metal-poor regions as bright (in ionizing photon density) as 
metal-rich region have, neverthless, 8 /im LSDs that are almost an order of magnitude lower. 
The most metal-poor regions can be as much as a factor of 30 fainter at 8 /xm than their 
metal-rich counterparts at fixed Paa LSD. In addition to the decrease of the dust opacity 
with metallicity (the same effect present for the 24 /xm), the data require a second ingredient 
to account for the underluminosity at 8 /im. The second ingredient, from the models of 
Draine &: Lil (120061 ). is the decrease of the mass fraction of the low-mass polycyclic aromatic 
hydrocarbon molecules in the dust mixture for decreasing metallicity. This decrease affects 
selectively more the 8 fim emission than the 24 /im emission, as the former has a larger 
fraction of its flux contributed by single-photon-excitation of PAHs. The models, indeed, 
account reasonably well for the observed deficiency of the 8 /im flux in metal-poor objects 
(IDraine et al.ll2007l ). A possible mechanism to reduce the mass fraction of PAH molecules in 
the low-metallicity galaxies is destr uction by the hard radiation field in those galaxies (e.g.. 
Madden et al.ll2006 : Wu et al. 2006 ). although these molecules have proven surprisingly ro- 
bust ( Engelbracht et alT2006 ). An alternative possibility to destruction is that the lowest 
metall icity galaxies may not have formed the carriers of the aromatic features in the first 
place jEngelbracht et all [2OO5I : bwekl 120041 ). 



The correlation that exists between LogfSs ^m, dust] and LoglS Paa,corr] for the metal-rich 
regions (equation 3) is statistically as significant as those existing for the 24 /im and for 
the combined optical/mid-IR luminosity. However, there are many caveats in using such 
correlation to trace SFRs with the 8 fim emission. First and foremost, stellar populations 
other than those which are currently forming stars contribute to the mid-infrared emission: 
thus, a calibration of the 8 /im emission as a SFR tracer will depend on the galactic area 
probed, and, possibly, on the ratio between current and past (or recent-past) star formation. 
A second limitation is the extreme sensitivity of the 8 /im emission to metallicity, which is 
about an order of magnitude larger than what observed for the 24 /im emission. 

In all cases, the presence of Active Galactic Nuclei in unresolved galaxies will complicate 
the use of the mid-IR band emission for SFR determinations. Galaxies in the SINGS sample 
that contain central non-thermal sources (Sy2s or LINERs) show prominent nuclear emission 
in the mid-IR bands; in our case, the 8 /im and 24 /im emission from the non-thermal nuclei 
tends to be as bright as some of the brightest HII knots in their host galaxies (Figures [3] and 

m). 
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9. Summary and Conclusions 

The analysis of the mid-infrared emission, at the Spitzer's IRAC 8 /im and MIPS 24 /im 
wavelengths, from a set of local galaxies drawn from the SINGS sample has shown that their 
viability as SFR indicators is subject to a number of caveats. The calibrations given in this 
work, provided the caveats are taken into account, should however be useful for measurements 
of SFRs in actively star-forming galaxies using their observed mid-infrared emission, in the 
absence of prominent central AGNs. 



The most robust of the indicators analyzed in this study is the one proposed by lKennicutt et al 



( l2007al ). which combines the observed Ha and 24 /im luminosities as probes of the total num- 
ber of ionizing photons present in a region. Here we present a calibration (equation 7) for that 
relation that is based on measurements of 220 HII regions/complexes in 33 nearby galaxies. 
Comparisons with models suggest that the calibration should be applicable to large systems 
or galaxies whose energy output is dominated by young stellar populations, up to at least 
luminosity surface densities of Spaa,corj- >10^°-10'^^ erg s~^ kpc~^. 

Second-best is the 24 /zm emission alone, which, although robust (within the limits of our 
analysis) for metal-rich objects, shows substantial deviations, at the level of factors 2-4, from 
the mean trend for decreasing metallicity. Models can account for the observed deviations 
as an effect of the increased transparency of the medium for lower metal abundances (less 
dust means less infrared emission). The calibration presented here (equations 8 and 9) is 
non-linear and reflects both observations and expectations from models. The calibration is 
appropriate for metal-rich HII regions or star bursts, and caution should be e xercised when 



applying it to extended star-forming galaxies (iKennicutt &: Moustakas II2006I ) or to metal- 
poor targets. In particular, galaxies with a prominent contribution from non-star-forming 
populations to the integrated light may receive a yet unquantified contribution to the 24 fim 
emission from a diffuse component. This diffuse component would not be directly related 
to star formation and may be in the single-photon-heating regime; its presence would add 
to the integrated 24 /im emission and would alter the calibration of this emission for SFR 
measurements. 

We do not present a calibration of the SFR based on the 8 /^m emission, as emission 
at this wavelength shows strong dependence not only on metallicity (more than an order of 
magnitude), but also on the size (about a factor of 2) and, possibly, star formation history 
of the region being measured. It should be remarked, however, that when measurements 
are limited to regions of star formation with a narrow spread in metallicity around the solar 
value, the 8 /xm emission shows a correlation that is almost linear with the ionizing photon 
rate. 
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The mean trends between the mid-IR emission (at 8 fim, 24 fim, and combined optical/mid- 
IR) and the ionizing photon tracer Paa are well accounted for by either a 4 Myr old burst 
of star formation or a >100 Myr old constant star formation model; this provides a measure 
of the degeneracies in the population models, but also a measure of the general applicability 
of the calibrations discussed here to both star-forming regions within galaxies and whole 
star-formation-dominated galaxies. 

In all cases, the spread around the mean trends has a r.m.s. of ~0.3 dex (each side). This 
spread is well accounted for by allowing stellar populations with a range of ages, between 
Myr and 8 Myr, to produce the observed emission. The older stellar populations can 
account for HII knots with high mid-IR emission relative to the ionizing photons density, 
although more complex dust geometries than those analyzed here can also contribute to the 
spread. 

Our analysis has concentrated on HII knots and starburst galaxies, where star-forming 
stellar populations dominate over more evolved populations in the bolometric output of 
the galaxy. However, the contribution of non-ionizing stellar populations to the heating of 
the dust emitting in the mid-IR region needs to be fully quantified, in order to test the 
applicability of (or derive modifications for) the SFR calibrations presented in this work to 
more general environments, such as quiescently star forming galaxies. This is the subject of 
a future investigation which era ploys the SINGS and other star-forming galaxies as testbeds 



flKennicutt fc Moustakas Il2006f ). 
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A. Models of Dust Absorption and Emission 

A.l. The Energy Balance 

To aid the interpretation of the observed correlations between the mid-infrared emission 
and the Paa emission (a tracer of the number of ionizing photons), we build a simple model 
of dust attenuation of the spectral energy distribution of a young stellar populations at UV- 
to-nearlR wavelengths, with the prescription that the dust absorbed energy is re-emitted 



in the mid/farlR with a spectral energy distribution that follows the iDraine fc Lil (120061 ) 
models. Ideally, the 'perfect' estimator of SFR should correlate linearly with the SFR itself. 
We have seen in section 6 that this is not the case for S24 and ^m, dust', a legitimate 
question is whether the nature of the observed non-linearity is such to jeopardize a reliable 
use of the mid-IR emission as a SFR indicator. Thus, our main interest in this section is to 
understand whether a heuristic model can account for the observed trends as a function of 
both luminosity surface density and metallicity, and is able to: 



• reproduce the observed mid-IR emission level and its trend as a function of the Paa 
LSD with simple prescriptions for the stellar emission spectrum, dust geometry, and 
ISM metallicity in the regions; 

• provide a physical explanation for the non-linear correlations between the mid-IR 
LSDs and the line emission LSD; 

• match overall expectations on the age of the stellar populations in the 'HII knots'. 



For the simplest ass umption that the dust extinction is described as an effective fore- 

ground attenuation A( A) (jCalzetti. Kinney fc Storchi-Bergmannlll994l ; lMeurer. Heckman fc Calzetti 
19991 ). the bolometric (3-1000 /xm) infrared luminosity L//j is given by: 



'IR 



POO 

/ F5(A)[l-10-°-^^(^)]dA, 
Jo 



(Al) 



where F5(A) is the stellar light SED. We use the 'foreground attenuation' as a working hy- 
pothesis, since it enables us to account for both absorption and scattering (in and out of the 

line of sight) of the stellar light by dust with a simple expression ( ICalzetti. Kinney fc Storchi-Bergmann 
I994J ; ICalzetti et al.ll2000l ). Potential limitations to this assumption will be discussed in sec- 
tion A.4. 

Our general results indicate a correlation between the mid-infrared emission and the 
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20061 ). and are qualitatively not dissimilar from the correlation between the bolometric 
infrared emission and the SFR derived for galaxies dominated by young stellar popula- 



tions using IRAS data (ILonsdale Persson fc Heloulll987l : lRowan-Robinson &: Crawfordlll989 



Devereux &: Younglll990l : ISauvage fc Thuanlll992l ). Equation Al shows that such correlations 



are mediated by the dust attenuation A(A); indeed, starburst and star-forming galaxies have 



crease ( 


Wang & Heckman 


1996 




Seckman et al. 


Moustakas, Kennicutt & Tremont 


i I2OO6 


). This g 



m our sample. A plot of the color excess b(B— VjU , derived from the tia/raa ratios of 
individual HII knots, as a function of the Paa LSD is shown in Figure [12], together with the 
best fit through the datapoints. The data in Figure [T2l are for the high metallicity subsample, 
the group with the largest number of datapoints and the widest range in color excess in our 
sample. If we introduce the ionizing photon rate per unit area, Ejon, in units of s~^ kpc~^ 
as derived directly from Spaa,corr for case B recombination, then in solar metallicity regions 
E(B— V) is related to the ionizing photon rate via the best fit (bi-linear on the log of the 
quantities) : 



E{B -V) = 0.21 



lO^^s ^kpc ^ 



{0.61±0.03) 



(A2) 



where E(B— V) is in magnitudes. The correlation has 5.7 a significance, and is independent 
of the specifics of the stellar population (age, burst or constant star formation, etc.). 

Some physical insights on the meaning of equation A2 can be gained by looking at the im- 
plica t ions of the scaling laws of star formation. The Schmidt Law (IKennicutt. Edgar fc Hodge 
19891 : lKennicuttlll998al ) provides a way to relate the SFR density to the gas density in galax- 
ies. Furthermore, in our own Galaxy, there is a tight correlation between th e color excess 

By combin- 



1978 



E(B— V) and the gas (HI+H2) column density (iBohlin. Savage Sz Drake 
ing the Schmidt Law, as derived locally for regions in NGC5194 (IKennicutt et al.ll2007al ). 
and the extinction-gas column density correlation, we get a relation between the color excess 
and the SFR density: 

E{B -V) = 10.30 7 J:^sf% (A3) 

where T^sfr is in units of Mq yr^^ kpc^^ and the factor 7 accounts for the fact that not all the 
gas (and the dust) is in front of the stars. Equation A3 can be readily related to equation A2, 
since the two exponents, 0.61 in the first and 0.64 in the second equation are formally the 
same number within the 1 a uncertainty of our fitting procedure (see, also. Figure [T^ . By 
fitting the data of Figure [T2] with a straight line of fixed slope 0.64 (equation A3), we get 
an estimate of 7 =0.4-0.45 for the constant star formation models, depending on the time 



^A(A)=k(A) E(B— V), where k(A) is the extinction curve, section 3.2 
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elapsed since the onset of star formation (6-100 Myr). This value of 7 is very close to 
the mean value of 0.5 expected in the case the star formation is located on average at the 
mid-point of the gas column density. 

Metallicity variations in the stellar populations produce small variations in equation A2. 
However, observations of the extinction to gas column density correlatio n for the Large and 



Smal l Magellanic Clouds s uggest roughly linear scaling with metallicity (IKoornneef fc Code 



198ll : iBouchet et al.lll985l ). This is likely to be the dominant dependence of equation A2 on 
metallicity, and we model it as a linear dependence on Z, the interstellar medium's metallicity 
in solar units: 

g(fl-V-) = 0.21z( ^^„J-^^_, ) . (A4) 

Although non-parametric (Spearman and Kendall) rank tests indicate that the data in 
Figure [12] deviate from the null (uncorrelated) hypothesis by 5.7 a, the spread around the 
mean trend is still significant. The 90% boundary is located at (-0.4 dex, +0.5 dex) along the 
vertical axis. We will see in section A. 4 that the large spread in Figure [T^ has a small impact 
on the main conclusions of this Appendix. That said, it should be important to keep in mind 
that although the trends have a physical base in the combination of the Schmidt Law with 
the extinction-gas column density correlation (Heckman, T.M. 2004, private communication; 
Tremonti, C.A., 2006, private communication), individual variations remain important. 

Equations A2 and A4 will be used in combination with equation Al in the following 
sections to provide a heuristic description of the relation between ionizing photon rates 
(star formation rates) and monochromatic dust emission in our regions, in order to un- 
derstand the broad trends observed in the data. Thus, equations A2 and A4 are derived 
from data and will be used to explain data, and the potential of circular argument should 
be addressed. The data in Figure [12] are completely independent of the data presented 
in Figures [3l H] and [7] The color excess E(B— V) is derived from ionized gas emission, 
while the monochromatic dust emission is mainly due to stellar continuum emission. Equa- 
tion Al also requires independent assumptions on the underlying stellar population, and 
the dus t geometry relative to both the stellar pop ulation and the ionized gas (see next 
section. ICalzetti. Kinney fc Storchi-Bergmannlll994l ). Therefore, the use of the extinction- 



ionizing photons correlation of equations A2 and A4 is unlikely to automatically enforce 
a fit to the data in Figures [31 [H and [7] Finally, we will see in section A. 4 that equa- 
tions A2 and A4 are not necessary to account for the observed data for Paa LSDs above 
~3xl0^^ erg s"^ kpc~^. 
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A. 2. Model Parameters 



For F.sfA) in equation Al, we adopt the stellar population models of Starburst99 (2005 



update; iLeitherer et al.lll999l ). Since we are exploring the infrared emission of both Hll 



knots and star-formation-dominated galaxies, we consider both bursts of star formation 
and constant star formation. Burst models are explored in the age range 0-10 Myr, as older 
ages do not produce enough ionizing photons to provide Paa-detectable HII knots in our 
data. As bursts age, the number of ionizing photons produced by the massive, short-lived 
stars, decreases, while the UV-optical stellar luminosity (the main contributor to the infrared 
emission) remains comparatively high, since it is contributed by lower-mass, longer-living 
stars. Thus, aging stellar populations produce high luminosity ratios L(IR)/L(Paa). For 
constant star formation, the representative case of a 100 Myr-duration model is considered. 
The difference between this model and a much younger case is the increase of optical and 
infrared light in the older model, due to the accumulation of previous-generations, low-mass, 
long-lived stars. The light from these older stars contributes to L(IR), but not to the ionizing 
photons budget. The net effect is that the 100 Myr constant star formation model produces 
about twice as much infrared emission as a, e.g., 6 Myr old case, for the same dust opacity. 
The variation between the 100 Myr and 1 Gyr models is about 6%. For the stellar IMF, 
we adopt the Starburst99 default, which consists of two power laws, with slope —1.3 in the 
range 0.1-0.5 Mq and slope —2.3 in the range 0.5-120 Mq. The SFRs derived from this IMF 
are a factor 1.59 smaller than those derived from a Salpeter IMF in the range 0.1-100 Mq, 
for the same number of ionizing photons. 

For the wavelength dependence of the dust a ttenuation in e quation A l, A(A)=k(A) 
E(B— V), we adopt the starburst obscurat i on curve ofICa ' 



1994; Meurer. Heckman fc Calzetti 1999; Calzetti et al 



zettil (I2OOII ) (see, also, ICalzetti. Kinney &: Storchi- 



2OOOI), which prescribes an attenua- 



tion of the stellar continuum a factor 0.44 that of the emission lines. 

The fraction of the infrared luminosity that emerges in t he IRAC 8 iim and MIPS 
24 /im bands, LfSl/Ld R) and L(24j / L(IR ), is from the model of IPraine fcH J2006h . which 
updates the model of iLi fc Draind (I2OOII ). At solar metallicity, the two papers provide 
consistent prescriptions for our case. The fraction of infrared light emitted in either 8 fim 
or 24 nm band inc reases non-linearly as a function of the starlight intensity (Figure 15 of 
Draine fc Lil (120061 )): the ratios are fiat at the low intensity end (single photon heating of 
dust) and rise sharply at higher intensities, producing 'hotter' infrared SEDs, i.e., peaked at 
shorter wavelengths. At solar metallicity, the dependence of L(8)/L(IR) as a function of the 
starlight intensity is much shallower than that of L(24)/L(IR); the former varies by 50% over 
the full range of in tensities analyzed h ere, while the latter changes by more than a factor of 
10. The results of iDraine fc Lil (120061 ) are expressed as a function of the parameter U, the 
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starlight intensity in units of the solar neighborhood value. 

We relate the typical stellar radiation field in our HII knots to U using the individual 
HII regions detected in the HST Paa image of the nearby (and, thus, better resolved), 
actively star-forming galaxy NGC6946. From the image we measure the sizes, defined as 
the radius of the circular area containing 80% of the Paa emission in each region, and 
line emission values of the HII regions. Total luminous energies are derived by associating 
the range of stellar models described above to the measured line intensities, and energy 
densities are derived by dividing for the measured volumes. By extrapolating the observed 
quantities, we find that burst models with mass 10^ Mq and ages <6 Myr or constant star 
formation models with SFR=1 Mq yr~^ (J-'Paa,corr ^10"^° erg s~^), with radii RsilOO pc, 
produce average starlight intensities U~400-1200. The factor ~3 variation in U accounts 
for two uncertainties in o ur derivation: (1) star forming reg ions have hotter SEDs than the 



local interstellar medium (IMathis. Mezger fc Panagialll983l ). and (2) our stellar models and 



regions' sizes cover a range of values/characteristics. The presence of dust absorption within 
the HII region and the conversion of some of the energy to free-free emission also contribute 
to the uncertainty (and some decrease) in U. However, the results described in the next 
section are fairly insensitive to the actual value of U, within the range of 3-4 uncertainty 
detailed above. The relation between number of ionizing photons and starlight intensity 
derived above suggests that, for constant radii of the HII regions, the transition from single- 
photon heating to thermal equilibrium heating for the 24 /xm emission begins to occur at 
emission line intensities Spaa,corr ~10^^ erg s^^ kpc~^. 

The impact of metallicity variations is explored by considering how the infrared emission 
changes under the two conditions that the regions/galaxies have solar metallicity {Z=Zq) 
and one-tenth solar metallicity (Z=l/10 Zq). The latter value is at the low end of the 
metallicity range in our sample; our goal is, indeed, to model the lower envelope to the data. 
The metallicity enters in two components of equation Al. The first component is the stellar 
population model, for which we use two metallicity values: solar and 1/5 solar (Starburst99 
does not provide models for our default 1/10 Z©, and 1/5 Z© is the closest metallicity 
value for the stellar populations). The second component is the dependence of the color 
excess on metallicity, which we describe in a linear fashion via equation A4. Additional 
metallicity-dependent ingredients are the conversion from ionizing photon rate to hydrogen 
nebular line luminosity (at the ~15% level; see section 3.2), and the frac tion of infrared light 



emerging in the 8 /xm and 24 /xm bands. We use, again, the model of iDraine fc Lil (120061 ) 
for the latter, with the assumption that variable metallicities have a major impact on the 
mass fraction of the polycyclic aromatic hydrocarbons in the dust mixtur^; this, in turn. 



Based on lDraine et al.l (|2007r) . a factor 9-10 difference in the mass fraction of small PAH is a reasonable 
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affects ttie dust emission in tlie mid-infrared bands, mainly tlie IRAC 8 /im band, wliicli 
results predominantly from single-photon-heating at low values of the starlight intensity. 
For reference, in the low-U regime, L(24)/L(IR) changes by about 60% and L(8)/L(IR) 
changes by an order of magnitude for a factor 10 var iation in the abuii dance of small PAH; 
the differences become negligible at the high-U end (jPraine fc Lill2006l ). 



A. 3. General Model Trends 

In the extreme case that we have no dust in our region, equation Al will be null, and 
there will be no relation between SFR (or ionizing photon rate) and dust emission. In the 
presence of even small amounts of dust, however, the relation established by equation Al will 
produce a general trend, with details that depend on a host of parameter assumptions. At 
the lowest levels of dust extinction, the relation between L(IR) and the number of ionizing 
photons is non-linear, with an asymptotic exponent of 1.64 for our equations A1-A4; a 
linear relation is established only once the integral on the right-hand-side of equation Al 
corresponds to most of the stellar energy. In other words, once the region contains enough 
dust that most of the stellar energy is absorbed and re-emitted in the IR, larger numbers of 
ionizing photons will linearly correspond to larger IR luminosities. 

The impact of the stellar population parameters on the characteristics of the 8 /im 
and 24 /xm emission has already been shown in Figures [9]-[TTl Stellar populations undergoing 
constant star formation since 100 Myr in the range SFR=4xlO~^-4 Mq yr~^ have equivalent 
characteristics to bursts of star formation of constant age 4 Myr and mass in the range 
M=10^-10^ Mq. For the burst models, the impact of age variation is investigated at constant 
cluster mass (10^ Mq) and in the two cases of constant or age-variable (equation A2) color 
excess. The general effect of increasing the age is to move the model curves towards smaller 
values of the ionizing photon rates for roughly constant IR emission, almost independently 
of assumptions on the variation (or constancy) of the color excess. 

The introduction of the age-variable extinction (via equation A2), however, attempts 
to mimic the observation that aging HII regions tend to be less dust extincted than the 
younger HII regions, whether because the more evolved populati ons have shed the nativ e 



cocoon or have migrated away from it through secular motions (IMayya fc Prabhul Il996l ). 
This also introduces a modest dependence of the IR SED and the fraction of infrared light 
emerging in the 8 /im and 24 fim bands on age, via changes in the starlight intensity that 
aging HII regions produce. However, the model does not include additional effects, such as 



assumption for metallicities between Z~0.1 and Z>0.5. 
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the expansion of HII regions as they evolve, that will affect the IR SED via the decrease of 
the dust temperature. 

Not surprisingly, metallicity variations have a strong effect on the observed mid-IR emis- 
sion as a function of the number of ionizing photons (Figure [T3|l . For the 24 /im emission, 
most of the effect comes from the fact that the ISM is more transparent at lower metallici- 
ties. For the 8 /im emission, roughly equal contributions are given by the more transparent 
medium and by the decreased mass fraction of low-mass PAH molecules in the dust mixture 
(Figure [131 left). 



The fiducial model for the IR SED as a function of starlight intensity (from lDraine &: Li 



20061 ) is compared with the more extreme assumption that the IR SED is constant for all 
analyzed regions/galaxies. For the assumption L(24)/L(IR)=const=0.3 (Figure \T3\ right), 
the model line has a different curvature from our fiducial one, starting with higher Log(S24^m) 
values at low Pact LSD, overshooting the datapoints in the LSD range Log(Spaa,corr)=39-40, 
and converging to a slope of unity at higher LSD values. Our fiducial model appears to 
better reproduce the slope of ~ 1.23 of the best fit to the datapoints (in the log-log diagram), 
and this slope is entirely due to the relation between L(24)/L(IR) and the starlight intensity. 
For the L(8)/L(IR) ratio, the change between a constant IR SED and a starlight-intensity- 
dependent SED is small, as the ratio changes by about 50% in the full range of LSD under 
consideration in this work. 



A. 4. Model Uncertainties and Limitations 

The simple assumptions on stellar populations and dust geometry made in the previous 
section lead to a number of limitations. Furthermore, the data show in some cases (e.g.. 
Figure [T2l) significant dispersion around the mean trend. The impact of considering different 
assumptions or including the full range (90th percentile) of data dispersion is briefly discussed 
here. 

The 90th percentile region for the E(B— V) versus Paa LSD correlation is shown in 
Figure[T2]by dotted lines that enclose 90% of the data around the mean correlation. The same 
region is reproduced on the 24 /im versus Paa LSDs in Figure [H] (left). We do not show the 
case for the 8 /xm LSD because the results are similar for both mid-IR bands. The dispersion 
in the color excess for fixed Paa LSD has minor impact on our baseline conclusions for the 
mid-IR emission, with significant impact (factor ~6 peak-to-peak variation) concentrated 
towards the low-intensity, and low extinction, regions. This effect is readily understood 
by recalling that once the color excess E(B— V)>0.5 mag, over 90% of the UV emission is 
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converted by dust into IR emission. Thus, the dependence of E(B— V) on the Paa LSD for 
Spaa,corr ^10^^ erg s~^ kpc~^ has httle impact on the amount of IR emission produced for 
stellar SEDs that emit most of their power in the UV. 

The independence of the IR luminosity on the color excess for high Paa LSDs is also 
shown for the case where E(B— V) is independent of Spaa,corr (i-e., equations A2-A4 are not 
applied). As a simplified example of E(B— V) independent of Spaa,corr, Figure [14] (right) 
shows the case of constant value E(B— V) = l for all Paa LSDs. As long as the color excess is 
sufficiently large that most of the UV stellar light is attenuated by dust, the model's trend 
for high luminosity regions does not change. However, when extended to low LSDs, the 
model with constant E{B—V) tends to exceed the low metallicity datapoints, and to show 
a flatter slope than the best flt line through the high-metallicity data. The latter is due to 
the models having reached the single photon heating for the 24 /xm emission, where there is 
a linear correlation between S24 and Spaa,corr- 

Age variations among different regions have been discussed in section 7, and they ef- 
fectively account for the spread around the mean trend for the infrared-versus-Paa LSDs 
(Figures IS HTT]) . Here we discuss the effect of adopting different dust models than our default 
one. 

Most non-foreground dust geometries, e.g., mixtures of dust and gas/stars, will have the 
effect of 'hiding' from direct detection a proportionally larger fraction of the stellar emission 
than foreground geometries, thus increasing the observed infrared emission, and S24^m and 
Ss/im, dust, relative to the ionizing photon rate {Spaa,corr) that can be recovered with simple ex- 
tinction correction methods (e.g., from Ha/Paa). Figure [H] (right) shows the effect of replac- 
ing our foreground dust assumption with an homogeneously mixed dust/stars geometry. The 
new model follows closely our baseline model up to LSDs Spaa,corr ~3xl0^^ erg s~^ kpc~^, 
and deviates upward of the baseline model for higher LSDs. This is the regime where the 
dust opacity is such that our standard procedure does not recover completely the number 
of ionizing photons and Spaa,corr is underestimated. If we were to allow for dusty cores with 
arbitrarily high attenuation values in our HII knots, the mixed model would deviate from 
the baseline model at arbitrarily low Paa LSDs, thus at least partially accounting for the 
dispersion of the datapoints around the mean trend. 

The competition of the dust with the gas for absorption of the ionizing photons will 
become increasingly more effective as the density o f the star formi n g regi ons increases, as in 



the case of LIRGs and brighter infrared galaxies. iRigby fc Riekd (120041 ) proposes that the 
lack of high excitation infrared flne structure lines in infrared galaxies is due to much of the 
massive star fo rmation occur r ing iri ultracompact HII regions, where such competition would 



be signiflcant. iDopita et al.l (120061 ) have modeled the behavior of such regions, conflrming 
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that the infrared continuum should get warmer, but also show the potential to suppress the 
emission lines substantially. 

One of the main assumptions in our analysis is that the 'typical' HIT knot in our sample is 
described by a single mean stellar population. Realistically, variations are expected, not only 
from galaxy to galaxy, but also from region to region. Under most circumstanc es, each of our 



apert ures will include multiple stellar populations covering a range of ages (ICalzetti et al. 



20051 ). Furthermore, the dust covering factor may change from population to population 
within each region, with older stellar populations often located in areas of lower extinction, 
or located further away from the dust they heat. The ratio of the 24 fim LSD to the observed 
Ha LSD, which by construction exacerbates any effect of dust extinction (Figure [T5l) can 
provide insights into this effect. The ratio 24 fim/Haobs spans more than two orders of 
magnitude for our data, with values between ~0.04 and ~10. The S24i_im/SHa,obs ratio is 
correlated (with a 7 cr significance for a Spearman non-parametric test) with the Paa LSD, 
with slope 0.83 ±0.03 in a log-log plane. This is a different way of expressing the correlation 
between star formation and dust extinction already shown in Figure | 



The model lines for solar metallicity are within the general locus of the data-points 
at low values of the Paa LSD, especially when the 90-percentile region from Figure [12] is 
included (Figure [T5l) . The brightest regions in Paa LSD, and most of the LIRGs, tend to 
be more consistent with a model where stars and dust are homogeneously mixed and there 
is no differential extinction between gas and star s (as assumed in all our baseline models. 



from lCalzetti. Kinney &: Storchi-Bergmamj|l994j ). This is not unexpected, as strongly star 



forming regions are i n general more heavily enshrouded in dust than less active regions 



( iGoldader et al.ll2002l ). In addition, in this regime, small errors in the measurement of the 
(faint) Ha line can produce large deviations of the data. The moderate discrepancy between 
the baseline models and data at the high luminosity end seen in Figure dS] has negligible 
impact on our results for the mid-IR LSDs as a function of the Paa LSD (section 6): in 
this regime, over 80% of the UV light is re-processed by dust into the infrared, and large 
fluctuations on the stellar light observed directly (including the nebular hydrogen lines) 
produce only small changes of L(IR). 

Despite the potential shortcomings of our assumptions, the ability of the models to 
describe the overall trends of the data enables us to dissect the individual 'ingredients' that 
produce the observed correlations to better understand their origin. 
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Table 1: Characteristics of the Sample Galaxies. 



Name 


Morph.'' 


Nucleus'' 


(km s^i) 


E(B-V)g" 


Dist.*^ 
(Mpc) 




12+log{0/H)'^ 


# Regions'^ 


HST ID' 


High Metallicity Galaxies 


NGC0925 


SAB{s)d 




553 


0.081 


9.12 


-20.33 


8.24-8.78 


10 


7919, 9360 


NGC1512 


SB(r)ab 




898 


0.011 


10.5 


-19.90 


8.37-9.05 


3 


9360 


NGC2403 


SAB{s)cd 




131 


0.043 


3.5 


-19.68 


8.31-8.81 


8 


7919 


NGC2841 


SA(r)b 


Syl/LIN 


638 


0.017 


9.8 


-21.12 


8.52-9.19 


5 


7919 


NGC2976 


SAc 




3 


0.074 


3.5 


-17.97 


8.30-8.98 


8 


7919 


NGC3184 


SAB{rs)cd 




592 


0.018 


11.10 


-20.46 


8.48-9.14 


6 


9360 


NGC3198 


SB(rs)c 




663 


0.013 


13.68 


-20.90 


8.32-8.87 


9 


9360 


NGC3351 


SB(r)b 




778 


0.030 


10.1 


—20.48 


8.60-9.22 


2 


9360 


NGC3627 


SAB{s)b 


Sy2/LIN 


727 


0.035 


8.7 


-21.17 


8.49-9.10 


6 


7919 


NGC3938 


SA{s)c 




809 


0.023 


12.2 


-20.10 


8.35-9.07 


9 


9360 


NGC4125 


E6 


LIN 


1356 


0.020 


21.4 


-21.89 


8.58-9.21 


1 


9360 


NGC4559 


SAB{rs)cd 




816 


0.019 


11.1 


-20.84 


8.25-8.79 


6 


7919 


NGC4569 


SAB{rs)ab 


Sy/LIN 


-235 


0.049 


16.58 


-21.90 


8.56-9.19 


5 


9360 


NGC4625 


SAB{rs)m 




609 


0.019 


9.17 


-17.53 


8.27-9.04 


6 


9360 


NGC4736 


(R)SA{r)ab 


Sy2 


308 


0.019 


5.3 


-20.59 


8.31-9.01 


5 


9360 


NGC4826 


{R)SA{rs)ab 


Sy2 


408 


0.044 


5.6 


-20.63 


8.59-9.24 


8 


9360 


NGC5033 


SA(s)c 


Syl.9 


875 


0.012 


13.28 


-20.87 


8.27-8.90 


9 


9360 


NGC5055 


SA(rs)bc 


LIN 


504 


0.019 


7.82 


-21.08 


8.42-9.13 


9 


9360 


NGC5194 


SA{s)bc 


Sy2.5 


463 


0.037 


8.2 


-21.43 


8.54-9.18 


43 


7237 


NGC5195 


SBO.l 


LIN 


465 


0.038 


8.2 


-19.99 


8.28-8.83 


1 


9360 


NGC5866 


S0.3 


LIN 


672 


0.014 


12.1 


-20.52 


8.43-9.02 


3 


9360 


NGC6946 


SAB{rs)cd 




48 


0.365 


5.0 


-21.11 


8.40-9.04 


9 


9360 


NGC7331 


SA(s)b 


LIN 


816 


0.097 


15.1 


-22.14 


8.40-9.05 


8 


9360 


Medium Metallicity Galaxies 


NGC1705 


SAO- 




632 


0.009 


5.1 


-16.29 


8.20-8.43 


3 


7919 


IC2574 


SAB(s)m 




57 


0.039 


2.8 


-17.33 


7.94-8.26 


6 


9360 


NGC4236 


SB(s)dm 







0.015 


4.45 


-19.12 


8.07-8.56 


3 


9360 


IC4710 


SB(s)m 




739 


0.095 


7.8 


-17.96 


8.11-8.62 


7 


7919 


NGC6822 


IB(s)m 




-57 


0.253 


0.47 


-14.97 


8.04-8.67 


3 


7919 


Low Metallicity Galaxies 


Ho H 


Im 




142 


0.035 


3.5 


-17.25 


7.68-8.07 


3 


9360 


DDO053 


Im 




20 


0.040 


3.56 


-13.68 


7.77-8.13 


6 


9360 


Ho IX 


Im 




46 


0.085 


3.3 


-13.67 


7.61-7.98 


3 


9360 


M81DwB 


Im 




347 


0.085 


6.5 


-14.20 


7.85-8.20 


4 


9360 


NGC5408 


IB(s)m 




506 


0.074 


4.8 


-17.22 


7.81-8.23 


3 


9360 


Discarded Galaxies 


NGC0024 


SA(s)c 




554 


0.021 


8.8 


-19.01 


8.32-8.92 




9360 


NGC1291 


(R.l)SB(l)0/a 




839 


0.014 


10.5 


-21.69 


8.42-9.01 




9360 


MSlDwA 


I ? 




113 


0.022 


3.55 




7.34-7.64 




9360 


NGC3034 


10 




203 


0.170 


5.2 


-20.51 


8.36-9.09 




7919 


NGC4631 


SB(s)d 




606 


0.018 


8.5 


-21.59 


8.13-8.76 




9360 


DD0154 


IB(s)m 




374 


0.010 


4.3 


-14.73 


7.54-8.02 




9360 



"Galaxy morphologies, nuclear activity, heliocentric velocity, and foreground Galactic color excess are from 
the NASA/IPAC Extragalactic Database (NED). The nuclear activity is reported for those cases where non- 
therr nal emission dominates the emission in the central region. The Galactic color excesses, E(B— V)g, are 
from Schlegel. Finkbeiner fc Davii ( 1998 ). 

* Adopted distances, in Mpc, as derived by Masters ( 2005[) . Note that changes in the galaxy distances do not 
affect the analysis, which is based upon luminosity surface densities (luminosity/area). 

'^Galaxy's absolute magnitude, based on Vj. from the RC3 (jDe Vaucouleurs et al.lll99ll ). as retrieved from 
NED. For NGC6822 the and for MSlDwB the m^ values have been used. 



'^Oxygen abundances; the two columns of values are from iMoustakas et al.l (|2007t ). see section 2 for a brief 
description of their derivation. 

'^Number of independent regions measured in the galaxy. 

^HST Program ID for the NICMOS observations used in this paper. 



Table 2: Galaxies Measurements. 



Name 


log Spaa,corr^ 


log ^Ha.obs 


[NHJ/Ha'' 


Av^ 


log Sgjim, duat'^ 


log S24Mm'' 




(erg s~^ kpc~^) 


(erg s ^ kpc ^) 


(mag) 


(erg kpc~^) 


(erg kpc~^) 


High Metallicity Galaxies 


NGC0925 


38.44±0.18 


39.18±0.08 


0.24 


0.45±0.19 


41.02±0.11 


40.48±0.11 


NGC1512 


38.55±0.06 


38.80±0.02 


0.47 


1.94±0.06 


41.09±0.06 


40.75±0.06 


NGC2403 


38.73±0.15 


39.27±0.08 


0.62 


1.07±0.17 


41.24±0.11 


40.35±0.11 


NGC2841 


37.83±0.30 


38.03±0.30 


0.55A 


2.12±0.43 


40.43±0.08 


40.08±0.08 


NGC2976 


38.50±0.15 


39.12±0.08 


0.30 


0.83±0.17 


41.12±0.10 


40.59±0.10 


NGC3184 


38.59±0.18 


38.72±0.11 


0.44 


2.31±0.21 


40.97±0.11 


40.69±0.11 


NGC3198 


38.56±0.18 


38.42±0.08 


0.36 


3.16±0.19 


41.20±0.10 


41.09±0.10 


NGC3351 


39.24±0.06 


39.48±0.04 


0.40 


1.97±0.07 


41.72±0.06 


41.68±0.06 


NGC3627 


38.92±0.23 


39.18±0.04 


0.55A 


1.94±0.23 


41.72dz0.11 


41.33±0.11 


NGC3938 


37.99±0.30 


38.76±0.06 


0.54 


0.38±0.31 


41.08d=0.08 


40.46±0.08 


NGC4125'^ 


37.13U 


36.00U 






40.15±0.06 


39.95±0.06 


NGC4559 


38.51±0.30 


39.33±0.08 


0.32 


0.24±0.31 


41.30±0.08 


40.75±0.08 


NGC4569 


38.29±0.30 


38.84±0.04 


0.50A 


1.05±0.30 


41.56±0.08 


41.35±0.08 


NGC4625 


38.06±0.11 


38.95±0.08 


0.46 


0.03±0.14 


40.79±0.08 


40.39±0.07 


NGC4736 


37.96±0.18 


38.88±0.04 


0.50<i 


0.01±0.18 


41.96±0.08 


41.54±0.08 


NGC4826 


39.28±0.06 


39.53±0.02 


0.55<^ 


1.97±0.06 


42.12±0.06 


41.69±0.06 


NGC5033'= 


38.57±0.08 ■ ■ ■ 






41.88±0.06 


41.35±0.06 




NGC5055 


38.61±0.11 


39.28±0.08 


0.50<^ 


0.68±0.14 


41.74±0.06 


41.17±0.06 


NGC5194 


38.95±0.06 


39.26±0.04 


0.50 


1.78±0.07 


41.54±0.06 


41.19±0.06 


NGCSigS'^ 


36.89U 


37.83U 






42.39±0.06 


41.35±0.06 


NGC5866 


37.64±0.30 


38.19±0.11 


0.55A 


1.04±0.32 


40.04±0.06 


39.47±0.06 


NGC6946 


39.44±0.06 


39.77±0.04 


0.56 


1.72±0.07 


42.34±0.06 


42.20±0.06 


NGC7331 


39.08±0.11 


39.14±0.08 


0.32 


2.54±0.14 


41.69±0.08 


41.20±0.08 


Medium Metallicity Galaxies 


NGC1705 


38.56±0.11 


39.48±0.11 


0.09 


0.03±0.16 


40.09±0.06 


40.06±0.06 


IC2574 


37.55±0.30 


38.44±0.08 


0.12 


0.14±0.31 


39.17±0.18 


39.18±0.18 


NGC4236 


37.57±0.30 


38.15±0.08 


0.17* 


1.09±0.31 


39.76±0.18 


39.42±0.18 


IC4710 


37.81±0.30 


38.56±0.08 


... g 


0.57±0.31 


39.64±0.11 


39.10±0.11 


NGC6822 


37.85±0.30 


38.74±0.11 


... g 


0.12±0.32 


39.46±0.30 


39.13±0.30 


Low Metallicity Galaxies 


Ho H 


38.37±0.30 


39.31±0.04 


O.lOA 


0.01±0.30 


40.01±0.13 


40.06±0.13 


DDO053 


37.97±0.30 


38.91±0.11 


0.04 


0.01±0.32 


39.33±0.18 


39.82±0.18 


Ho IX^ 


36.12±0.73 


36.32±0.2 


0.05A 




36.92U 


38.14±0.30 


MSlDwB 


37.76±0.18 


38.60±0.08 


0.04 


0.29±0.19 


39.03±0.11 


39.15±0.11 


NGC5408 


38.26±0.30 


39.20±0.04 


0.02* 


0.01±0.30 


39.29±0.06 


40.05±0.06 



"Average luminosity surface density and extinction in the central ~50"x50" of each galaxy (~144"xl44" for 
NGC 5194). The extinction-corrected and the 'observed' (not extinction corrected) values are listed for Paa 
and Ha, respectively. 

''[NII](A 6584 A)/Ha, as obtained from spectroscopy ( Moustakas et al.ll2007 ). and before convolution with 
the narrowband filters' transmission curves. An 'A' after a value indicates an adopted (non measured) 
value, for those cases where a spectrum is not available or the available spectrum is dominated by a central 
non-thermal source. Adopted values come from galaxies of comparable metallicity to the target ones. 

'^For NGC4125 and NGC5195, the only detected sources are the central Sy2 nuclei; for these, only 8 /im and 
24 /zm emission is detected, while Ha and Paa are upper limits. 

'^For NGC4736, NGC4826, and NGC5055, the [NH]/Ha values are derived from the comparison of the HST 
and ground-based narrowband images. The HST narrowband filters centered on Ha reject almost completely 
[NH], thus providing a reference for the ground-based images. 

'^For NGC5033, no Ha data are available. 

•^V alues from the spectrum of IHo. Filippenko k SargentI UMI) (NGC4236) and 



Masegosa. Moles fc Campos-Aguilail (jl9M) (NGC5408). 



^For IC4710 and NGG6822, no ratios are available. However the [NH] contribution to the narrowband Ha 
filter is negligible for these two galaxies, since both galaxies have low metallicity and the narrowband filters 
transmit less than 3% and 4% of the light from the 6548 A and 6584 A [NH] lines, respectively. 

''The Paa emission for HoIX is the middle value between the upper limit measured from the HST/NICMOS 



- 52 - 



Table 3: Characteristics of the Starburst Galaxies. 



Name 


Morph.^ 




E(B-V)g" 


Rb 


12+log{0/H)'= 






(km s~-^) 




(arcsec) 




UGCA292 


ImlV-V 


308 


0.016 


5.1 


7.2 


SBS0335-052 


BCG 


4043 


0.047 


4.1 


7.3 


HS0822+3542 


BCG 


732 


0.047 


4.1 


7.4 


VIIZw403 


Pec. 


-103 


0.036 


6.1 


7.7 


UM461 


BCD/Irr 


1039 


0.018 


5.1 


7.8 


Mrkl450 


Comp. 


946 


0.011 


5.1 


8.0 


IIZW40 


BCD 


789 


0.820 


11.2 


8.1 


NGC5253 


Im pec 


407 


0.056 


14.2 


8.2 


NGC2537 


SB(s)m 


431 


0.054 


8.1 


8.7 


NGC2146 


SB(a)ab 


893 


0.096 


13.8 


8.4-9.0 



° Galaxy morphology, heliocentric velocity, and foreground Galactic color exce ss are from the NASA/IPAC 



-V)g, 



is from Schlegel. Finkbeiner fc David 



Extra galactic Database (NED). The Galactic color excess, E(B 
(|l998h . 

''Radius, in arcseconds, of the region of active star formation, as measured in the NICMOS images. 

'^Oxygen abundance s, reproduced fr om Table 1 of lEngelbracht et al. For UGCA292, the oxygen 

abundance is from [PLl vugin ("20011), for NGC5253 from iMartiiJ ilMZl); for NGC2146 the line ratios of 
Ho. Filippenko fc Sargent (.IQQTi) have been converted to a range of possible oxygen abundances using the 



strong lines method of iKewlev et al.l (|2002| ) 



Table 4: Measurements of the Starburst Galaxies. 



Name 


log Spaa^ 

(erg kpc~^) 


log Sgfim, dust'^ 

(erg kpc~^) 


log S24Aim'' 

(erg s"-"^ kpc~^) 


UGCA292 


38.44 


40.11±0.33*' 


39.71±0.42 


SBS0335-052 


39.59 


41.67±0.02 


41.97±0.10 


HS0822+3542 


39.12 


39.48 ±0.21 


40.68±0.11 


VIIZw403 


39.2 


40.48±0.24 


41.34±0.11 


UM461 


39.42 


40.57±0.09 


41.43 ±0.10 


Mrkl450 


39.72 


40.73±0.09 


41.64±0.10 


IIZw40 


40.28 


41.59±0.05 


42.4±0.10 


NGC5253 


40.37 


42.34±0.07 


42.93±0.10 


NGC2537 


39.94 


42.02±0.16 


41.94±0.10 


NGC2146 


40.64 


43.44±0.03 


43.24±0.10 



"The luminosity surface density at Paa, 8 /im, and 24 /im in the starburst regions, with radius listed in 
Table [3l The Paa LSDs are from the HST SNAP program 9360, and are only corrected for foreground 
Galactic extinction (Table [3]). Measurement uncertainties for Spaa are around 15%-20%. The LSDs in 
the Spitzer bands are 'whole galaxy' measurements, corrected to infinite aperture, although in all cases the 
central starburst (measured in the NICMOS images) is the dominant contributor to the flux. 

''The 8 una emission from this galaxy is affected by a latent image from a previous observation. Every effort 
has been made to remove the contaminating latent image from the measurement, but the presence of some 
small remnant contamination cannot be excluded. 
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Fig. 1. — Example of aperture selection in one of our fields. The central ~1' of the galaxy 
NGC0925 is shown at Ha (top-left), Paa (top-right), stellar-continuum-subtracted 8 fim 
(bottom-left), and 24 /im (bottom-right). The 13" apertures used for photometric measure- 
ments are shown as red circles, and correspond to physical sizes of ~580 pc. The field shown 
is one of the two obtained in the central region of this galaxy. North is up; East is left. 
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Fig. 2. — Histogram of the multiplicative factors applied to the observed Paa line emission 
to correct for the effects of dust extinction, as determined using the Ha/Paa ratio (sec- 
tion 3.2). The vertical axis shows the number of HII knots to which each correction factor 
is applied. The high (12+log(0/H)>8.35), medium (8.00<12+log(O/H)<8.35), and low 
(12+log(O/H)<8.00) metallicity regions (section 2) are shown as three separate histograms, 
coded as red dot-dash line, green continuous line, and blue dash hne, respectively. Most 
corrections are less than 50%. 
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Fig. 3. — (a) Luminosity surface density at 8 /im, Sg ^m, dust, as a function of the extinction- 
corrected LSD at Paa, Spaa,corr, for the 220 HII knots in the 33 galaxies for which photometric 
measurements have been obtained. LSDs are averaged over 13" photometric apertures. 
The 8 fim emission is stellar continuum-subtracted (section 3.1). Data points are divided 
into three metallicity bins: high (red filled triangles), medium (green stars), and low (blue 
asterisks) oxy gen abundance (sec t ion 2) . Filled black squares mark the local starbursts from 
the sample of lEngelbracht et al.l (120051) (section 5.1). 3 a error bars are indicated for each 
data point. Open black star symbols indicate the location of the non-thermal sources (Sy2s 
or LlNERs, section 4.1) and open black circles indicate extended background sources. The 
best fit line through the high metallicity (red) datapoints is shown as a continuous line, 
while the dashed line is the linear fit through the same datapoints with fixed slope of 1. (b) 
Histogram of the deviation of the HII knot data in panel (a) from the best fit line through 
the high metallicity data (the continuous line in panel (a)). The values of the best-fit 
coefficients are c=(0.94±0.02) and d=(4.80±0.85) (equation 2). Three separate histograms 
are shown, for high (red), medium (green), and low (blue) metallicity data. The medium 
and low metallicity histograms have been multiplied by a factor 2 to make them visible. 
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Fig. 4. — (a) and (b). The same as Figure [3l for the luminosity surface density at 24 /im, 
S24 ^lm■ In addition to the same datapoints as Figure [3l panel (a) als o repo rts the Luminous 
InfraRed Galaxies (LIRGs) from the sample of lAlonso-Herrero et al.l (120061) (black asterisks; 
section 5.2). The values of the parameters (c, d) in the horizontal axis of of panel (b), are 
given in equation 3, and are c=(1.23±0.03) and d=(— 6.88±0.97). 
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Fig. 5. — The ratio of the mid-IR to the extinction-corrected Paa LSDs, as a function of 
the galaxies' distances, for the 8 /xm (left) and the 24 /im (right) HII knot measurements, 
respectively. Background sources and non-thermal sources are not included. In both panels, 
the high-metallicity datapoints are consistent with no correlation of the ratio as a function of 
distance, implying that the correlations between the mid-IR emission and the Paa emission 
are not driven by distance effects. The 8 /im emission of the medium and low-metallicity 
HII knots remains on average deficient relative to that of the high-metallicity data also when 
only galaxies at similar distances are considered (thus flux measurements are performed in 
similar-size regions). With the possible exception of NGC 6822, which is at a distance of 
only 0.47 Mpc and shows lower-than-average values for its metallicity bin, the observed 
8 fim emission deficiency in metal-poor regions is not an effect of a bias in the size of the 
regions that are being measured. 
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Fig. 6. — Ratio of the 8 /im to 24 /xm LSD as a function of the extinction-corrected Paa LSD 
for the 220 HII knots and for the local starbursts. Symbols and colors are as in Figure [31 
The deficiency of the 8 /xm emission, relative to the 24 /xm emission for the low-metallicity 
datapoints is independent of the average ionizing photon rate in the region. The decrease of 
the 8 /zm to 24 /im LSD ratio as a function of increasing Paa LSD for the high metallicity 
points indicates that the component of thermal equilibrium dust contributing to the 24 /xm 
emission is increasing in stren gth (the dust is in thermal equ ilibrium and 'warmer' at higher 
ionizing photon densities, see, lHeloulll98d : braine fc Ljbood ) . A contribution to the decrease 
of the 8 /im emission due to increased dest ruction rate of the carr iers for increasing starlight 
intensity (Paa LSD) may also be present (IBoulanger et al.lll988l ). 
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Fig. 7. — (a) and (b). As Figure HI for the linear combination of Haobs and 24 fim LSD 
(equation 4). Symbols are as in Figures [3] and HI (Panel a): Data for the HIT knots and 
the LIRGs. The best fit line through the high metallicity HII knots (continuous line) is not 
significantly different from a linear relation with slope of unity (dashed line). (Panel b): 
Histogram of the deviation of the HII knot data in panel (a) from the best fit line through 
the high metallicity data (the continuous line in panel (a)). Unlike Figures [3] and HI the 
histograms of the medium and low metallicity datapoints have not been multiplied by a 
factor 2. The values of the x-label parameters {c, d) are derived from equation 4 and are 
c=(1.02±0.02) and d=(-0.74±0.97). 
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Log[Sp„ ^^^J (erg/s/kpc2) 

Fig. 8. — Mid-IR or combined-optical/mi d-IR LSD as a function of the extinction-corrected 
Paa LSD, averaged over the entire central region covered by the HST/NICMOS observations, 
for each of the star-forming galaxies from Table 1 (section 4.1). The three panels show on 
the vertical axis, from top to bottom, the 8 /xm, 24 /im, and the linear combination of Ha 
and 24 /xm LSDs in logarithmic scale, as indicated by the label at the bottom-right corner 
of each plot; the vertical axis has the same units as the horizontal axis, erg s^^ kpc~^. Each 
HST/NICMOS image is about 50" in size, except for NGC5194, where a region of 144" in 
size has been observed. Color coding of each galaxy is the same as the HII knots in Figure [31 
The continuous lines are the best linear fit through the high metallicity (red triangles) data, 
after excluding the Sy2-dominated fluxes of NGC4569, NGC4736, and NGC5195 (marked 
as black stars), thus leaving 19 independent datapoints. The dotted lines are the best fits 
through the high metallicity HII knots from Figures [3l HI and [7] for Ss^m, dust, S24^m, and 
a Sjj„. „(,„+b So/i,.™. resDectivelv. 




Log[Sp„ ,„,J (erg/s/kpc2) 

Fig. 9. — Luminosity surface density at 24 /im as a function of the extinction-corrected 
Paa LSD for the same datapoints as Figure H] (after removal of the Sy 2 nuclei, the fore- 
ground and background sources, and the NGC5033 HII knots data, leaving 164 independent 
datapoints in the high metallicity HII knot subsample). The continuous line shows the 
best linear fit through the high metallicity HII knots, from Figure HI Models of infrared 
and ionized gas emission are superimposed on the data, for a variety of star formation 
histories, stellar population ages, and metallicity (see Appendix). Models with solar metal- 
licity {Z=Zq) ism and stellar populations include: 100 Myr-old constant star formation 
(SFR/area=4xlO~^-4 Mq yr~^ kpc~^, long-dash line); instantaneous burst with variable 
mass (lO'^-lO* Mq kpc~^) and color-excess, and constant age of 4 Myr (dot-dashed line); 
instantaneous bursts with constant mass (10® Mq kpc~^) and variable age, and both variable 
color excess (dot-dashed line with empty circles) and constant color excess (E(B— V)=2 mag, 
dot-dashed line with empty squares). The circle and square symbols mark the population 
ages, right-to-left: 0.01, 2, 4, 6, 8, 10 Myr. The dotted line marks a 1/10 Zq model of 
constant star formation over the past 100 Myr. The upward-pointing arrow marks the 
approximate luminosity where the transition between single-photon heating and thermal 

pniiiliKriiiTn hpatincr for thp Hiist Vipcrins t.n nppiir 
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Fig. 10. — Luminosity surface density at 8 /im as a function of the extinction-corrected Paa 
LSD for the same regions/galaxies as Figure [9] (minus the LIRGs). The continuous hne is 
the best hnear fit through the high metaUicity HII knots, as in Figure [31 Models of infrared 
and ionized gas emission are the same as Figure M 
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Fig. 11. — Combined-optical/mi d-IR LSD as a function of the extinction-corrected Paa 
LSD for the same data as Figure [91 The continuous hne is the best hnear fit through the 
high metalhcity HII knots, as in Figure [3 Models of infrared and ionized gas emission are 
the same as Figure O 
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Fig. 12. — The color excess E(B— V), measured from the Ha/Paa ratio (section 3.2), as a 
function of the Paa LSD, for the 164 HII knots in the high metalhcity subsample. The best 
bi-hnear fit through the data (continuous hne) has slope 0.61 (equation A2). The dotted 
lines mark the region containing 90% of the datapoints around the best fit line. The dash 
line is from equati on A3, where the slope is kept at the fixed value 0.64 (from the Schmidt 
Law in NGC5194. [Kennicutt et aPboOTah . 
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Fig. 13. — Models of dust absorption and emission for constant star formation populations 
and variable metallicity are compared with data, for the 8 /xm (left) and the 24 /xm (right) 
LSD as a function of the Paa LSD. Symbols for the HII knots, the local starbursts, and the 
LlRGs are as in Figures [914T01 Straight continuous lines are the best fits through the high- 
metallicity HII knots, similar to the lines in Figures [3] and |H All models are for a 100 Myr old 
constant star formation population, and include: solar metallicity (Z=Z0, long-dash line); 
1/10 Zq and standard low-mass PAH molecules fract ion (short-dash lin e); 1/10 Z© and 
depleted low-mass PAH molecules fraction (dotted line lDraine fc Lill2006l ). For the 24 /im- 
versus-Paa plot, the effect of variations in the IR SED are also explored; in particular, 
our default assump tion for the shape of the IR SED as a function of the starlight intensity 
( iDraine fc Lill2006l ) is compared with the extreme assumption that the IR SED is constant, 
i.e., L(24)/L(IR)=const=0.3 (dot-dash line). 
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Fig. 14. — The 24 /xm LSD as a function of Paa together with the best fit hne (continuous 
straight hne, Figure HI). Left. The fiducial model of a dusty 4 Myr old stellar cluster 
with increasing mass (dot-dashed line) is bracketed by the dispersion curves derived from 
including in equation Al the 90 percentile region of the E(B— V)-versus-Spaa,corr correlation 
(dotted lines and Figure [T2|) . Right Model lines for a 4 Myr old stellar cluster of increasing 
mass, and two assumptions for the dust distribution that are different from our baseline 
model: (1) foreground dust geometry and constant E(B— V) = l as a function of Spaa,corr 
(short-dash line), and (2) a homogeneously mixed dust-star geometry with no differential 
extinction between gas and stars, and variable E(B— V) according equation A2 (long-dash 
line). 




Fig. 15. — The ratio of the 24 /xm to the observed Ha LSD as a function of Paa, for the 
HII knots and LIRGs (Figure [9]), compared with models. The datapoints are correlated with 
a ~7 a significance; the continuous line is the best bi-linear fit through the data, and has slope 
0.83. Solar metallicity models include (Figure [9]): instantaneous burst with variable mass 
(10^-10^ Mq) and color-excess, and constant age of 4 Myr (dot-dashed line); instantaneous 
burst with constant mass (10^ Mq), variable age and variable color excess (dot-dashed line 
with empty circles marking, right-to-left: 0.01, 2, 4, 6, 8, 10 Myr). The dotted lines mark 
the upper and lower boundaries to the 90-percentile region from Figure [121 The long-dash 
line is the model of homogeneous dust-stars mixture, with no differential extinction between 
gas and stars, from Figure [TH The short-dash line is the 1/lOth solar metallicity model, 
and marks to lower envelope to the datapoints. The upward-pointing arrow marks the 
approximate luminosity where the transition between single-photon heating and thermal 
equilibrium heating for the dust begins to occur (Figure [9l). 



